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ELECTRON  GAS  BEHAVIOR  IN  THE  BASE  AND  MIXING 
REGIONS  OE  A  WEAKLY  IONIZED  PLASMA  FLOW  AT  AN 
ABRUPT  INCREASE  IN  CROSS  SECTION 

By 
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March,  1968 

Chairman:   Dr.  R,  B.  Gaither 

Major  Department:   Mechanical  Engineering 

An  experimental  investigation  is  made  of  the  electron 
behavior  in  a  weakly  ionized  argon  plasma  jet  (with  an 
initial  boundary  layer)  exhausting  into  a  quiescent  medium. 
The  low  pressure  (1  to  20  mm  Hg)  subsonic  jet,  studied 
under  conditions  v/aich  are  steady  but  not  necessarily 
field  free,  is  found  amenable  to  hoth  theoretical  and 
experimental  analyses. 

The  flow  field  of  the  plasma  jet  is  artificially  broken 
down  into  three  separate  regions:   the  potential  core 
region,  the  jet  mixing  region,  and  the  developed  flow 
region.   Results  obtained  in  the  core  and  developed  flow 
regions  are  in  agreement  with  those  obtained  by  other 
investigators.   Particular  attention  is  given  to  the 
electron  gas  behavior  in  the  base  and  jet  mixing  regions 
due  to  both  normal  and  superimposed  potentials  in  the  base 
region.   A  mathematical  model,  based  on  a  known  initial 
electron  density  profile,  is  developed  for  predicting  elec- 
tron number  densities  within  the  jet  flow. 
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Detailed  velocity,  electron  temperature,  and  electron 
number  density  profiles,  obtained  using  conventional  gas 
dynamic  probing  techniques,  including  the  Langmuir  probe, 
are  presented. 
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1.   INTRODUCTION 

In  recent  years  the  particular  result  of  an  electrical 
discharge,  the  creation  of  a  plasma,  has  received  consider- 
able attention  (_1)  (2^)  (3^).  Notable  among  the  engineering 
applications  spurring  on  these  recent  studies  are  plasma 
propulsion  systems,  rocket  exhaust  and  re-entry  simulators, 
magnetohydrodynamic  generators,  energy  conversion  devices, 
arc  welders,  etc. 

While  basic  research  in  plasmas  has  been  and  continues 
to  be  attractive  to  a  large  number  of  investigators  (4)  (5) 
(6^),  research  with  a  clear  objective  of  obtaining  information 
in  flow  systems  particularly  significant  to  engineering- 
design  remains  somewhat  limited.  This  is  particularly  true 
in  situations  where  the  fluid  dynamics  of  a  problem  are  not 
fully  understood. 

A  flow  situation  of  this  type,  which  is  of  considerable 
interest  in  gas  dynamics,  and  consequently  in  plasma  dynamics, 
is  found  in  pressure  induced  flow  systems  where  an  ionized 
gas  proceeds  past  an  abrupt  increase  in  cross  section  into 
a  zone  where  for  some  distance  the  flow  is  not  fully 
developed.   This  includes  flows  in  such  common  configurations 
as  free  jets,  ejector  systems,  and  flows  past  orifices.   In 
all  of  these,  the  behavior  of  a  neutral  fluid  in  the  base 
and  mixing  regions  is  both  difficult  to  predict  and  of 
considerable  engineering  significance.   If  the  fluid  is  a 
plasma,  the  situation  is  further  complicated. 


■'■Underlined  numbers  in  parentheses  in  the  text  of  this 
dissertation  refer  to  references  listed  in  the  b  i  hliogmiihy . 
Numbers  in  parentheses  refer  to  equations. 
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The  objective  ^..  ^..iv^:^  ^^.^.-  -i>.^_  jl- 

mentally  determine  the  ochavioi"  o: 
ponents  in  a  weakly  ionized  gas  proc 
velocities  past  an  abrupt  increase  in  crc. 
quiescent  re':::icn.   Both  the  f;-as 

ties  were  l.-. ,  ^.. ...  ^^„ 

correlating-  charged  particle  diffusion,  ui.  j^iy 

rates  to  the  macroscopic  >: 
attention  was  attached  to  studying  c 
behavioi'  in  the  base  and  jet 
by  an  applied  potential  ir 

Investigations  c-  ^..^  >,^.^  ^.j. ...... ^-^  ^^■^^.■^^^..^^. 

sor.ic  jet  flows  have  been  made  by  Kuethe  (7) 
and  Squire  (9^)  ,   Using  the  concept  of  an 
fluid,  they  measured  such  properties  as  p 
length,  velocity  decay,  quantity  of 
Their  experimental  findings  arrse: 
theoretical  results.   ^....^1.-:.  ^-  ^- 

similar  solutions  to  "che  incc 
equations,  one  subgroup  of  v/hiv_, 
symmetric  laminar  jet.   Howeve:- 
valid  only  far  from  the  noziile,  i.e. 
region. 

Very  few  v.'or.ii^  .- 
gas  dynamic  mixing  of  jvC   flow^ 
distances  fi-om  a  no.:zle  >:.xit. 
contributions  are  the  v. 
Kleinstein  (13) .  Chov 

developed  soiutior...  -..^     .-..^.._      ...^... 

laminar  coir.pressible         whil.      

primari: 

develop.  ohnin, 

flov/s      jrupt  increases  in  cros. 

channel. 
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In  the  more  complex  problem  of  a  plasma  flow  near  a 
nozzle  exit,  the  current  information  available  is  even 
more  restricted  than  it  is  in  the  pure  gas  dynamic  problem. 
Research  in  this  area  has  been  limited  to  a  great  extent 
by  the  lack  of  suitable  diagnostic  methods  and  the  prohibi- 
tively high  temperatures  involved.   Nevertheless,  several 
noteworthy  contributions  have  been  made  in  this  area.   Using 
Langmuir  probe  techniques,  Gaither  (15)  investigated  the 
electron  gas  behavior  in  the  potential  core  region  of  a  low 
pressure,  field-free  argon  jet  flow.   Graf  (16^)  used  both 
microwave  and  Langmuir  probe  techniques  to  investigate  the 
electron  density  distribution  in  the  fully  developed  region 
of  a  low-density,  free  expansion  argon  jet.   Other  signifi- 
cant studies  have  been  made  by  Warder  (17^)  ,  Sherman  (_18)  , 
and  Cobine  (19).   These  studies  will  be  considered  in  more 
detail  later  in  this  dissertation. 

In  the  following  chapter  some  of  the  more  fundamental 
aspects  of  plasma  behavior,  especially  those  of  importance 
in  the  configuration  under  study  here,  are  presented. 
These  basic  considerations  form  a  background  for  Chapter  3 
which  deals  with  the  theory  of  the  experiment.   Chapters  4 
and  5  are  concerned  with  the  experimental  procedure,  results, 
and  comparison  of  these  results  with  theory.   In  Chapter  6 
the  conclusions  drawn  from  this  experiment  are  presented, 
and  recommendations  for  future  investigations  are  set 
forth. 


2.   THEORETICAL  CONSIDER, 

2 . 1  Plasmas 

An  ionized  gas  is  essentially  i.      dus  r.. 

ions,  fi'ee  electrons,  and  neutral  ■>;   . 

stances  where  the  ionized  gas  is  e.cccrically 
significant  lengths  in  the  system  are  _ 
radii  of  the  effective  force  fields  surrc 
particles,  the  ionized  gas  is  often  called  a  p 

In  the  inert  and  noble  gases.  cc.;r.r.:onlv  r.~. 
studies,  negative  ions  cannot  I; 
electrons  to  the  neutral  molecu-.e. 
ment  and  charge  transfer  reactions  are  ause.. 
argon-^,  an  inert  gas,  is  used  : 
assumed  that  the  plasma  flows -bei:.        -ered 
negative  ions . 

It  is  conv^  .^..  „  ..  „  ^_.^o   _...         ...^  _ 

property  which  will  be  ;     ^ulai- 

discussions  u.\l.    ,/h        osse;L, 

within  the  gaseous  mixture.   Th     oparty  is  cal 

degree  of  ionization  (15)  and  is  defined  by 


ox  :\ 
■  Con. 


j9. 
foi- 


,  O  —  V^  . 
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where  n^  is  the  ion  number  density  and  n^  is  the  neutral 

particle  number  density.   This  relation,  as  expressed  here, 
is  restricted  in  utility  to  the  case  of  single  ionization 
and  for  neutral  plasmas  having  n^C:^   n^,  where  n^  is  the 
electron  number  density.   Since  it  is  possible  to  have 
plasmas  of  low  degrees  of  ionization  as  well  as  plasmas 
which  are  completely  ionized,  £  may  vary  from  zero  to  one. 
It  is  important  to  note  that  e  is  a  relative  term  and  does 
not  depend  upon  absolute  values  of  the  numerical  density 
terms.   Therefore,  it  is  possible  to  have  plasmas  of 
different  densities,  exhibiting  completely  different 
properties,  with  the  same  degree  of  ionization.   However, 
there  are  many  similarities  which  exist  between  low  and 
high  density  systems,  and  the  investigations  of  low  density 
systems  where  prohibitively  high  temperatures  can  be 
avoided  often  lead  to  generalized  results  which  can  be 
applied  to  high  density  systems  (21) . 

Plasmas  may  be  produced  through  energization  by  elec- 
tromagnetic, mechanical,  nuclear,  or  thermal  means. 
Depending  upon  the  method  by  which  the  ionization  is 
produced  and/or  maintained,  the  gas  may  be  in  either  an 
equilibrium  or  non-equilibrium  state.   In  the  equilibrium 
state,  the  temperature  of  the  gaseous  constituents  are 
approximately  equal  and  extremely  high. 

Non-equilibrium  plasmas  are  usually  produced  by  methods 
other  than  thermal  and  at  low  pressures  and  densities  where 
the  gaseous  constituents  can  be  maintained  at  different 
temperatures.   In  low  density,  electrically  energized 
plasmas,  the  electron  temperature  is  usually  of  the  order 
of  thousands  of  dcgroos,  while  the  ion  and  neutral  particle 
temperatures  are  kept  at  values  of  a  few  hundred  degrees. 
The  non-equilibrium  cliaracter  of  the  plasma  is  determined 
primarily  by  the  .strength  and  configuration  of  the  external 


field  producing  the  ionization,  as  well  as  by  the  energy 
transfer  mechanisms  between  the  constituents. 

In  a  DC  excited,  low  pressure  gas  system,  such  as  the 
one  under  consideration,  the  plasma  is  in  a  non-equilibrium 
state.   Only  the  charged  particles  receive  energy  directly 
from  the  electric  field.   The  electrons,  because  of  their 
low  mass  and  consequent  high  mobility,  receive  sufficient 
energy  to  impart  to  them  high  velocities  and  consequently 
high  temperatures.   The  ions,  on  the  other  hand,  because 
of  their  relatively  large  mass  are  not  as  greatly  affected 
by  the  applied  field,  at  least  in  a  direct  sense. 

One  very  convenient  method  of  experimentally  studying 
the  non-equilibrium  plasma  as  it  tends  toward  a  state  of 
thermal  equilibrium  is  to  use  an  electrode  configuration 
which  permits  the  plasma  to  flow  to  a  region  remote  to  that 
of  the  applied  external  field.   As  the  plasma  is  caused  to 
leave  the  discharge  region,  the  collisional  processes  in 
the  plasma  tend  to  cause  it  to  approach  a  state  of  thermo- 
dynamic equilibrium.   As  pointed  out  by  Goldstein  (21),  in 
the  absence  of  applied  fields  the  interaction  processes  that 
lead  to  this  equilibrium  state  are  the  case  of  all  the 
phenomena,  whether  observable  or  not,  that  take  place  in 
the  plasma. 


2.2   Jet  Flow  Structure 

The  electrode  configuration  chosen  for  this  investigation 
is  shown  schematically  in  Figure  1.   It  consists  simply  of  a 
stick  cathode  and  an  annular  anode.  This  particular  configura- 
tion has  the  inherent  advantage  of  forcing  all  the  gas  to 
flow  through  the  established  electric  field  before  proceeding 
into  the  test  section.   A  glow  discharge  is  created  in  the 
plasma  generator,  and  the  positive  column  of  this  discharge 


7 
.is  forced  through  the  anode  as  a  result  of  the  prossrro 
difference  established  between  the  generate^  -....  .„^.  _^c- 
tions . 

The  flow  field  of  the  resulting  axially  symmetric  jet 
may  be  artificially  broken  down  into  three  distinct  regions: 
the  potential  core  region      jet  mixing  region,  and  the 
developed  region.   These  regions,  shov/n  in  Fi'":v.re  2,  3. 
analyzed  independently  wi...  ...^  ^w....^..^.,  ^^..^^  ^..^..._.   .. 

matched  at  the  interface  of  the  assigned  regie, 
the  coupled  analysis  is  valid. 

A  jet  flow  leaving  the  nozzle  at  subsonic  velocities 
contains  a  potential  core  of  length  Lq,    composed  oJ  the  gas 
of  the  original  stream.   This  gradually  narrowing  subsonic 
core  is  suri-ounded  by  a  gradually  widening  jex   ..._.>.^-, 
region,  an  annular  region  in  which  the  original  gas  s  li-c..: 
and  the  surrounding  quiescent  fluid  have  intermixed. 
Mixing  continues  to  extend  across  the  jet  cross  section 
until  the  potential  core  disappears  at  length  L^ ;  the: 
the  flow  is  developed,  and  mixing  t.kes  place  throurrhout 
the  jet . 

In  the  po;;eii";;_aj.  co-e  region  of  xho      ^.ne  Luj^k  or 
macroscopic  properties  remain  relatively 
this  portiOxi  of  the  flow  field  may  Le  : 
as  inviscid.   The  general  conical  shapo  of  this  1 
experir-cntallv  vorifiod  bv  sevcrr^l  i. 

the  coro  and  the  qi.^  ^  flu;             .jro^ccu. 

axial  direction.  ore  radius  r^  h. 

decrease  linearly  with  axial  distanc 

that  it  may  be  matho...  Hy  descr:. 


/  .1 


^c  =  ^o<-^  -- 


L^ 
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where  Rq  is  the  nozzle  radius.   Clearly,  this  equation  is 
suitable  only  when  the  initial  boundary  layer  at  the  nozzle 
exit  is  negligibly  small,  and  the  jet  flow  leaves  the 
nozzle  at  a  uniform  velocity  Uq .   Otherwise,  R  is  replaced 
by  an  effective  nozzle  radius  R^ff,  over  which  the  axial 
velocity  is  found  to  be  uniform. 

The  axial  length  of  the  core  region  at  low  subsonic 
velocities  for  an  axially  symmetric  plasma  flow  was 
experimentally  verified  by  Gaither  (1_5)  to  be  between  4  and 
5  nozzle  diameters.   One  of  the  objectives  of  Gaither 's 
investigation  was  to  determine  the  core  length  and  its 
effect  upon  the  mixing  of  the  plasma  components  with  the 
surrounding  gases  with  attention  being  given  to  the  effect 
of  core  length  upon  plasma  properties. 

In  the  core  region  of  an  ionized  gas  flow,  a  very 
significant  effect  which  occurs  is  that  there  is  a  much 
greater  lateral  diffusion  by  charged  particles  than  by  the 
neutral  particles.   This  is  primarily  due  to  action  by  the 
electron  gas  component  which  possesses  an  extremely  high 
mobility  (due  to  light  mass  and  high  temperature) .   The 
ion  gas  component  follows  in  accordance  with  the  ambipolar 
effect  to  be  discussed  later.   Recombination  docs  not 
become  a  significant  factor  until  the  electron  temperature 
nears  that  of  the  ions  and  neutral  particles  (23^)  (2_4)  . 

At  the  point  where  the  emerging  jet  flow  and  the 
quiescent  fluid  meet,  a  free  jet  boundary  is  formed.   The 
zone  between  the  free  jet  boundary  and  the  potential  core 
is  called  the  mixing  region.   This  is  a  zone  of  complicated 
mass  and  energy  transfer,  and, at  a  sufficiently  high  Reynolds 
number,  it  becomes  turbulent  in  character.   The  effects  of 
this  mixing  region  on  the  behavior  of  the  charged  particle 
components  in  an  ionized  flow  had  not  been  determined  prior 
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to  this  investigation.   The  lack  of  such  information 
initiated  the  particular  study  described  herein. 

The  base  region,  as  shown  in  Figure  2,  is  that  region 
immediately  outside  the  mixing  region  and  near  the  nozzle 
exit . 


2.3   Ionization  and  Recombination 

Although  the  region  of  most  interest  in  this  investi- 
gation is  the  plasma  stream  after  it  leaves  the  discharge 
region,  it  is  necessary,  in  view  of  future  discussions,  to 
consider  briefly  the  ionization  processes  occurring  in  the 
plasma  generator.   The  consequent  recombination  processes 
will  also  be  discussed. 

The  energy  necessary  for  ionization  may  be  imparted  to 
an  atom  as  it  passes  through  a  discharge  by  an  electron 
collision,  a  positive  ion  collision,  a  collision  with  a 
high  energy  neutral  atom,  the  absorption  of  a  quantum  of 
radiant  energy,  or  by  some  combination  of  these.   However, 
at  gas  temperatures  equivalent  to  less  than  one  electron 
volt,  the  most  important  mechanism  for  increasing  the 
charged  particle  population  in  a  gas  subjected  to  an  applied 
electric  field  is  ionization  by  electron  impact.   A  neutral 
particle  will  be  ionized  by  an  approaching  electron  if  the 
combined  energies  of  the  two  particles,  relative  to  the 
motion  of  the  center  of  mass,  exceed  the  ionization  poten- 
tial of  the  gas  V^^ ,  i.e., 

■^rrj^lVe-Ql  +  -f-'^A,/^  -^  >  Vi^  (3) 

where  m^ ,  mj^,  v^,  and  v,^  are  the  electron  and  noulr.nl 
particle  masses  and  velocities,  respectively  (25)  (26).   G 
is  the  center  of  mass  velocity  defined  by 
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Elimination  of  G"  between  equation  (3)  and  (4)  yields  the 
following  criterion  for  ionization  in  a  single  collision: 

-jrMirk^J%-%l^>V^'  (5) 

Since  m^yy  xciq   and   j  v^  -  v^jt^i  v^,  this  criterion  reduces 
to  the  simple  relation 


Lvi^V^>V-'  (6) 

Thus,  the  electron  must  acquire  a  kinetic  energy  in  its 
acceleration  by  the  applied  electric  field  which  is  at  least 
as  great  as  the  ionization  potential  for  the  neutral  gas. 

The  first  ionization  potential  for  argon  Vi-.,  i.e., 
the  energy  required  to  remove  the  first  electron  from  an 
argon  atom,  is  15.8  ev,  while  the  second  ionization  poten- 
tial Vi^  is  27.6  ev  (27)  (28).   Each  succeeding  ionization 
potential  is  correspondingly  higher.   The  probability  of 
ionization  higher  than  first  order  in  systems  at  electron 
temperatures  such  as  the  one  under  consideration  is  very 
low.   Integration  of  the  Maxwellian  energy  distribution  from 
15.8  ev  to  oo  shows  that  less  than  0.2%  of  the  electrons 
possess  energies  sufficient  to  cause  ionization  for  an 
argon  plasma  maintained  at  3  mm  Hg  and  exhibiting  an 
electron  temperature  of  approximately  5,000  K.   Under  the 
same  conditions  only  0.005%  of  the  electrons  possess 
energies  greater  than  27.6  ev.   Thus,  it  is  clearly  seen 


that  only  single  ioniza . 
in  this  investigatxoix. 

As  shown  in  the  preceding  discussion,  very  slow  movii.:-' 
electrons  do  not  produce  ionization.   Very  high  energy 
electrons  are  also  relatively  poor  ionizers  since  thej"^ 
pass  through  an  atom's  sphere  of  influence  withou 

an  electron  (19).   Thus,  maximum  ionization  -.  ,.. ^^ 

within  an  intermediate  range  of  electron         :,  alt. 
ir  must  be  pointed  out  that  this  optimum      /  for  ioniza- 
tion greatly  exceeds  those  atta     o  in  this  investigation, 
In  argon^  average  electron  kinetic      ^ies  of  1  ev  (about 
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11,000  K)  give  rise  to  ionizing  collision  probabilities*^ 

of  4/cm.   A  maximum  is  noted  w- \:  ^.^^^■^.^:. -^-„,  — .«  -.,  -- 

approximately  1.35  x  10-^  K,  where  the  probability  ■=  ~ 
S3 /cm  (2S_)  . 

The  preceding  discussion  applied  only  to       -se  of 
ionization  by  single  collisions        :  pos; 
electron  having  -    ..ere:'''  less  than  that  re: 
ionization  to  excite  -  .-^..„-^;^  ^ — „_w..^  \.^   — .  -_.. 
level  above  che  ,;i"-ound  state,  i.e.,  to  ar.  - :--■ 
The  neutral  p:ir'oicle  may  the/.     -onized  by  =ubs«^ 
collisions      energetic  electrons.   This  : 
is  usually  referred  to  as  cumulative  ioni. 

A  large  br.      of  atons  ir.  circited  state:  is  ^^rc: 

very  short  '  '  the  order  of  10~^  sc_o-_i,. 

there  are  ce.        :ited  states,  usually 
ij.etas table  energy  .        where  the  probability  ^ 
transition  is  extrenely  small  or  is  prohibited  by 


'^Collision  p.^obability  as  us 
av:  ■   ■  colli 

truv.      ^  ^,..  ^..^'ourh  ^....  ^...\,v, ., 
Tn  -  . 
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selection  rules  of  quantum  mechanics.   Argon  has  two  such 
metastable  states  at  11.5  and  11.7  ev,  respectively  (18) . 
Metastable  states  occur  principally  in  noble  gases  and  have 
mean  lifetimes  of  the  order  of  milliseconds.   With  such 
comparatively  long  mean  lifetimes,  atoms  in  metastable 
states  may  undergo  several  collisions  before  radiating 
spontaneously.   Thus,  they  may  lose  their  excitation  energy 
through  collisions  without  radiating.   Such  collisions  are 
usually  i-ef erred  to  as  collisions  of  the  second  kind. 

According  to  recent  studies  by  Brewer  and  McGregor  (29) 
on  low  density,  arc-heated  argon  flows,  there  are  large 
concentrations  of  metastable  atoms  present  in  all  such 
plasma  flows.   However,  experimental  information  in  this 
area  is  lacking,  and  no  analytical  techniques  are  currently 
available  that  permit  accurate  determination  of  the  number 
densities  of  metastable  atoms  in  such  flow  systems.   The 
most  promising  technique  appears  to  be  through  the  use  of 
optical  absorption  measurements. 

In  the  flow  region  downstream  of  the  nozzle  only  weak 
field  forces  are  present.   This  is  due  to  the  particular 
electrode  configuration  being  used  and  the  fact  that  a 
plasma  is  incapable  of  supporting  a  large  potential  gradient 
(usually  less  than  1  volt/cm).   Thus,  further  ionization  is 
essentially  negligible.   Experimental  verification  of  this 
will  be  presented  later. 

Allowing  that  further  ionizing  effects  are  not  present 
in  this  region,  the  rate  of  change  of  ion  and  electron 
number  densities  is  simply  the  rate  of  recombination  given 
by 

dlt^  ^dUL.^  -^v^n-r::f<^v^i  (7) 
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where  the  quasineutrality  condition  n^  c^  n^    is  imposed. 
cK  is  the  coefficient  of  recombination.   It  is  a  function 
of  electron  temperature,  as  well  as  electron  density,  gas 
pressure,  and  type  of  gas  involved  (30). 

Experimental  values  of  c<  are  difficult  to  obtain 
directly  in  jet  flow  systems  because  the  electrons  are  lost 
from  the  jet  by  both  diffusion  and  recombination.   Conse- 
quently, most  measurements  of  o<  for  electron-ion  recombina- 
tion are  made  using  microwave  techniques  and  probes  in  the 
afterglow  of  non-flow  decaying  plasmas  (28)  .   The  reader 
is  referred  to  the  works  of  Biondi  (3^^)  and  Bates  and 
Dalgarno  (32)  for  recent  developments  in  experimental 
techniques  and  values  for  recombination  coefficients  in 
neon  and  argon.   Biondi  (31.)  observed  a  decaying  argon 
discharge  (P  in  the  range  of  1  to  25  mm  Hg ;  T^  C^   300  K; 
n   of  the  order  of  10^  to  lO^/cm^)  and  found  the  electron- 
ion  recombination  coefficient  to  be  approximately  10"   cm^/ 
sec . 

There  are  several  known  mechanisms  for  electron-ion 
recombination.   The  more  important  ones  are  enumerated  as 
follows:   (a)  radiative  recombination,  (b)  dielectronic 
recombination,  (c)  dissociative  recombination,  and  (d)  three 
body  recombination.   In  the  following  paragraphs  these 
mechanisms  are  considered  separately,  with  comments  as  to 
their  applicability  to  this  investigation. 

In  radiative  recombination,  a  free  electron  and  a 
positive  ion  undergo  a  collision  in  which  radiative  emis- 
sions of  energy  occur  as  given  by  the  following  reaction 
equation 


■h   e" >-  A  ■*-  ^^' 


(8) 
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h  is  Planck's  constant  and  i>  is  the  frequency  associated 
with  the  emitted  photon.   This  reaction  has  a  very  small 
cross  section  and  is  normally  of  no  significance,  except 
in  plasmas  of  very  low  densities  (18)  (28)  .   According  to 
Papoular  (33)  ,  the  total  effective  cross  section  (5"^  for 
such  recombination  decreases  as  \/vq    .      He  gives  a  table 
of  C^    values  for  several  low  electron  energies,  where  the 
radiative  recombination  coefficient  c/^  is  given  by 

oC^:^^^^-  (9) 

Considering  a  typical  set  of  conditions  for  the  discharge 
undei"  investigation  where  an  average  value  of  Tg  is  5,000  K, 
Vg  is  determined  to  be  1,2  x  10^  m/sec .   <5^  is  approximately 
1.5  X  10-19  cm3/sec  (33).   This  gives  a  value  of  1.8  x  10-19 
cm^/sec  for  c?<^  ,  which  is  negligibly  small  compared  to  the 
recombination  coefficient  of  10-6  cm^/sec  given  by  Biondi  (31) 
for  a  similar  argon  discharge.   Thus,  it  is  highly  improbable 
that  this  type  of  recombination  could  be  present  to  any 
significant  degree  in  the  plasma  flow  system  under  study. 

Dielectronic  recombination  is  represented  by  the 
reaction  equation 

h     -h    Q   ^A       .  (10) 

In  this  type  of  recombination  reaction,  an  electron  and  a 
positive  ion  undergo  a  collision  in  which  the  two  are 
bound,  with  part  of  the  excess  energy  being  absorbed  by  a 
second  electron  in  the  neutral  atom.   This  produces  a  doubly 
excited  neutral  atom  A**  which  is  unstable.   Under  circum- 
stances where  reionization  does  not  occur,  the  excited  atom 
then  proceeds  to  complete  the  recombination  process,  either 
by  undergoing  a  radiative  transition  to  a  lower  energy  level 
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or  by  transferring  the  excess  energy  to  other  particles 
during  the  collision  processes.   Dielectronic  recombination 
coefficients  are  believed  to  be  generally  smaller  than 
corresponding  values  for  radiative  recombination  (18)  (28) . 

Traces  of  rare  gas  molecules,  e.g.,  Aq ,  in  argon 
discharges  have  been  verified  by  mass-spectrographic  studies 
by  Tuxen  (34) .   The  possible  presence  of  such  molecules  in 
argon  plasma  flows  as  well  as  static  discharges  leads  to 
consideration  of  the  mechanism  of  dissociation  recombination. 
This  mechanism  can  be  described  for  argon  by  the  reaction 
equation 


i^'  (11) 


The  electron  initially  collides  with  the  molecular  positive 
ion  A2,  transforming  it  into  an  excited  state.   A  radiation- 
less  transition  then  occurs  in  which  the  molecular  bond  is 
broken,  and  the  atoms  of  the  excited  molecule  recede  one 
from  the  other.   Further  de-excitation  may  then  take  place 
(18)  (28).   Although  possibly  present,  dissociation  recom- 
bination has  not  been  reported  to  be  of  much  importance  in 
rare  gas  discharges  (18). 

Three  body  recombination  can  be  represented  by  the 
reaction  equation 

A  -^g'-^-X 5^ /4 -^X  .  (12) 

If  an  electron  encounters  a  positive  ion  in  the  vicinity  of 
a  third  body  X,  it  may  impart  sufficient  energy  to  the 
third  body  to  slow  it  down  to  where  it  will  readily  recom- 
bine.   The  third  body  may  be  another  particle  or  a  containing 
wall.   Tlie  possibility  also  exists  that  turbulence  in  the 
mixing  region  of  a  .jot  flow  could  provide  a  psuodo  third  body, 
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In  gas  discharges  near  and  below  10  mm  Hg»  it  has  been 
found  that  the  dominant  mechanism  for  removing  the  charged 
particles  from  the  plasma  is  ambipolar  diffusion  (to  be 
discussed  later)  (35) .   Neutralization  of  these  particles 
usually  takes  place  at  the  container  walls  by  three  body 
recombination.   Further  comments  regarding  possible  recom- 
bination mechanisms  in  the  plasma  flow  of  this  study  are 
presented  in  Chapter  5, 


2.4   Equilibrium  and  Continuum  Considerations 

The  plasma  model  adopted  in  this  investigation  is  that 
of  a  set  of  charged  gases  (electron  and  positive  ion)  and 
a  neutral  gas  mutually  penetrating  into  each  other.   Par- 
ticles of  each  gas  interact  with  the  other  particle  types 
and  with  each  other.   The  transfers  of  momentum  and  energy 
during  these  interactions  are  the  basic  microscopic  events 
governing  the  plasma's  macroscopic  properties  and  leading 
to  the  transport  phenomena  of  diffusion,  mobility,  viscosity, 
and  thermal  and  electrical  conduction. 

During  collisions  appreciable  kinetic  energy  exchanges 
are  made  only  between  particles  of  relatively  equal  mass. 
Thus,  in  a  plasma  the  mean  energy  of  the  electrons  and  the 
mean  energy  of  the  heavy  particles  are  only  weakly  coupled. 
Under  circumstances  where  the  velocities  of  all  of  the 
particles  in  a  plasma  are  allowed  time  to  randomize,  the 
heavy  particles  quickly  achieve  an  equilibrium  distribution 
with  temperature  Tj^  (or  Tj^ ,  assuming  Tj^  —  Tj^)  after  only 
a  few  collisions  with  each  other.   The  electrons  also  achieve 
an  equilibrium  distribution  with  temperature  T^  after  a  few 
electron-electron  collisions.   Thus,  for  some  period  of 
time  two  temperatures  exist.   In  the  absence  of  any  selective 
accelerating  field  they  will  approach  a  final  equilibrium 
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temperature.   The  number  of  collisions  required  for  Tg  and 
Tj-j  to  relax  to  an  equilibrium  temperature  is  approximately 
equal  to  the  neutral  particle  to  electron  mass  ratio,  i.e., 
approximately  mj^/mg  collisions  (2)  . 

Complete  thermodynamic  equilibrium  exists  in  a  plasma 
system  only  if  all  of  the  components  of  the  system  are  in 
thermal,  chemical,  electrical,  and  mechanical  equilibrium 
with  each  other  as  well  as  within  their  own  species.   It 
is  recognized  that  this  is  not  the  case  in  the  jet  flow  field 
under  consideration  since  finite  property  gradients  exist 
in  all  flow  systems.   However,  since  the  flow  is  subsonic 
and  the  particle  densities  are  high  enough  to  support 
continuum  analyses,  these  gradients  are  very  small.   The 
deviation  from  local  equilibrium  is  then  small  (15)  (18)  . 
The  effect  of  the  large  difference  between  the  electron  and 
neutral  particle  temperatures  will  be  shown  later  to  have 
negligible  effect  on  the  macroscopic  character  of  the  flow. 
Thus,  the  picture  of  the  plasma  being  studied  is  that  of  a 
mixture  of  gases  with  each  component  having  an  equilibrium 
(or  at  least  quasi-equilibrium)  distribution,  characterized 
by  a  temperature  T.   The  densities  are  such  that  close  energy 
coupling  between  the  neutral  and  ion  species  exists,  but 
the  coupling  between  the  electrons  and  heavy  particles  is 
weak. 

The  assumption  of  radiation  equilibrium  is  implicit 
in  the  assumption  of  local  thermodynamic  equilibrium. 
Radiation  originates  in  the  system  as  a  result  of  recombina- 
tion of  charged  particles  and  de-excitation  of  excited  atoms. 
However,  the  energy  contribution  due  to  radiation  in  systems 
at  extremely  low  degrees  of  ionization  is  negligibly  small 

(18)- 

As  a  basis  for  future  discussion,  it  is  important  to 

examine  the  component  ga.ses  of  a  plasma  flow  to  see  whether 
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they  behave  as  continuous  media  or  exhibit  characteristics 
having  to  do  with  molecular  structure.   The  Knudsen  number 
K  is  a  measure  of  the  importance  of  these  molecular  effects 
and  is  defined  by 

K=-^»  (13) 

where  7\     is  the  mean  free  path  of  the  component  gas  and 
L   is  some  characteristic  length  of  the  system.   A  Knudsen 
number  of  0.01  is  usually  taken  as  the  value  above  which 
molecular  effects  should  not  be  ignored  (36),  i.e., 

-A_<O.Oi  (14) 

Lt 

is  taken  as  justification  for  using  a  continuum  analysis 
to  describe  a  flow  situation. 

For  a  continuum  boundary  layer  theory  to  be  applicable 
in  the  jet  mixing  region  for  both  charged  and  neutral  gases, 
it  is  necessary  that  the  appropriate  mean  free  path  be 
small  compared  to  the  mixing  region  thickness  6    ,  or 

-ZL-<0.02.  (15) 

For  this  particular  investigation  the  values  of  both 
?\     and  S     vary  with  position.   Thus,  values  of  K  must  be 
examined  near  the  center  of  the  jet  flow  and  near  the  outer 
edge  of  the  jet  mixing  region. 

According  to  simple  kinetic  theory  (no  field  or  body 
forces  being  considered)  in  the  scattering  of  the  ith-type 
particle  by  the  jth-type  particle,  the  mean  free  path  is 
given  by 
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where  the  first  subscript  refers  to  the  particle  being 
scattered,  and  the  second  subscript  refers  to  the 
scatterer  (27^).   Qi_-j  is  the  cross  section  for  momentum 
transfer,  and  n^  is  the  number  density  of  the  jth-type 
particle.   The  1/2   does  not  appear  in  the  denominator  of 
equation  (16)  for  electron-heavy  particle  interactions 
because  of  the  large  mass  differences.   A  method  for 
determining  the  necessary  mean  free  paths  has  been  out- 
lined by  Graf  (16)  and  is  given  in  Appendix  B. 

Another  method  for  justifying  a  continuum  flow 
hypothesis  has  been  presented  by  Gaither  (15) .   Following 
a  procedure  outlined  by  Shapiro  (37)  and  Jeans  (38) ,  he 
showed  that  if  the  boundary  layer  thickness  S     was  the 
smallest  characteristic  dimension  of  the  system,  i.e., 
L  =  ^  ,    the  criterion  for  continuum  flow  was  given  by  the 
inequality 


^L> 


M 


100.  (17) 


Both  the  Reynold's  number  Rg  and  the  Mach  number  M  are  to  be 
determined  at  the  point  of  interest  in  the  flow  field. 


2.5  Plasma  Properties  and  Equation  of  State 

The  properties  used  to  describe  an  ionized  gas  system 
can  be  conveniently  classified  into  two  categories: 

1.  State  properties  which  describe  the  system  locally 

2.  Transport  properties  which  describe  the  behavior 
of  the  system  as  regards  diffusion,  thermal  and 
electrical  conduction,  and  viscosity. 
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2,5.1   State  Properties  and  Equation  of  State 

The  use  of  argon  as  a  test  gas  at  low  pressures  (1  to 
10  mm  Hg)  and  at  temperatures  highex*  than  room  temperatures 
(approximately  300  K)  permits  the  assumption  of  ideal  gas 
behavior  (27).   However,  justification  must  be  given  to 
permit  treatment  of  the  charged  particle  components  as 
ideal  gases  because  of  the  interparticle  forces  involved. 
Requiring  that  the  inertial  forces  be  much  greater  than 
the  Coulombic  forces,  Gaither  (15)  showed  that  the  charged 
particle  components  may  also  be  treated  as  ideal  gases  if 
the  following  criterion  is  satisfied: 


3     3 

Ti  <.<  7.3  T   X  10  . 


(18) 


T  and  n  are  the  temperature  (in  degrees  Kelvin)  and  number 
density,  respectively,  of  the  component  under  consideration. 
As  may  be  seen,  this  inequality  is  easily  satisfied  for 
experimental  conditions  where  the  number  densities  are  less 
than  10l2/cm*^  and  the  corresponding  temperatures  are  at 
least  a  few  hundred  degrees. 

Allowing  that  each  of  the  component  gases  can  be 
treated  as  an  ideal  gas,  the  corresponding  equations  of  state 
are  given  by 


■Pj=n3^T;, 


(19) 


where  j  denotes  the  jth  component,  and  k  is  Boltzmann's 
constant.   P.;  is  the  partial  pressure  of  the  component,  and 
Tj  denotes  the  translational  kinetic  temperature. 

Assuming  the  ideal  gas  relation  and  Dalton's  law  of 
partial  pressures  to  hold  for  the  type  of  gas  system  under 
consideration,  the  total  pressure  of  the  system  may  be 
expressed  as 
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■?=Zn'^^^^^>  (20) 


or 

P=  H-n^T^  +  -n-Ji  +  ri^T<z ).  (21) 

For  a  weakly  ionized  gas  at  low  pressure,  where  nj^»  n-|^ , 
Tj^CiTj^,  rigCii  n^,  and  T^  <  10^  K,  the  last  two  terms  in 
parentheses  on  the  right  hand  side  of  equation  (21)  are 
negligibly  small  compared  to  the  first  term.   Thus,  the 
equation  of  state  under  such  conditions  is  given  by 

P=-n^-feT,j.  (22) 

It  is  important  to  note  at  this  point  that  although 
the  gas  flow  is  slightly  ionized,  this  ionization  is  not 
sufficient  to  appreciably  affect  the  macroscopic  character 
of  the  flow.   Thus,  the  flow  may  be  adequately  described 
in  terms  of  the  state  properties  of  the  neutral  component. 

Other  state  properties  may  be  determined  by  applying 
the  appropriate  relations  for  mixtures  of  ideal  gases. 
These  properties  can  also  be  approached  with  the  knowledge 
that  the  charged  particle  components  have  very  little 
effect  on  their  behavior  if  the  degree  of  ionization  is 
very  small. 


2.5.2   Transport  Properties 

The  transport  properties  of  a  gaseous  medium  are  those 
properties  that  determine  the  behavior  of  the  gas  as  regards 
diffusion,  mobility,  thermal  and  electrical  conduction,  and 
viscosity.   Transport  properties  arise  from  the  fact  that 
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the  component  particles  of  all  fluid  media  are  in  motion 

and  constantly  interacting  with  each  other.   Unlike  state 
properties,  the  transport  properties  are  related  to  the 
momentum  transport  cross  sections  between  the  particles. 
These  cross  sections,  in  turn,  are  usually  complicated 
functions  of  temperature  and  the  interaction  mechanisms. 
The  current  difficulties  involved  in  the  determina- 
tion of  transport  properties  are  given  by  Sherman  (18)  as 

(a)  lack  of  a  suitable  theory,  (b)  lack  of  appropriate 
physical  section  data,  and  (c)  absence  of  techniques  for 
carrying  out  the  complex  calculations. 

Two  approaches  are  commonly  taken  in  determining 
transport  properties:   (a)  the  exact  or  formal  method,  and 

(b)  the  approximate  method.   The  so-called  "exact"  method 
starts  with  the  Boltzmann  equation,  and  the  mathematical 
difficulties  involved  are  enormous.   Thus,  only  under 
certain  assumptions  and  for  special  cases  can  any  useful 
solution  be  obtained  (2^)  .   For  the  purposes  of  this 
investigation  only  approximate  methods  will  be  considered. 

Chapman  and  Cowling  (39)  invoked  the  approximate  but 
simple  concept  of  the  mean  free  path  /\  in  establishing 
the  transport  coefficients  of  a  single  component  monatomic 
gas.   They  expressed  the  coefficients  of  diffusion,  viscos- 
city,  and  thermal  conductivity  as  follows: 

3)=-^'^'^^  (23) 


yU.  =4^1_>?777>l  V,  (24) 
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li=^f-V^^'^'>  (25) 
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where 

Vz 


In  these  equations  v  is  the  mean  particle  speed,  n  is  the 
number  density,  m  is  the  particle  mass,  k  is  Boltzmann's 
constant,  T  is  the  temperature,  and  Q  is  the  momentum 
cross  section. 

Fay  (2^)  suggested  that  the  above  equations  could  be 
extended  to  the  case  of  a  partially  ionized  gas  by 
assuming  that  the  transport  of  mass,  momentum,  and  energy 
within  the  plasma  could  be  expressed  as  the  summation  of 
the  transports  due  to  each  particle  species.   He  suggested 
the  following  summation  formulae: 

Dj  =  f^^i^i'  (28) 


|J.^^.^r,,'mi■^,^i^ 
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(29) 


*-   iZQ  L 


(30) 


where         \-   ,^  " -^rrr'  (31) 

J 


Q^_-j  is  the  momentum  transport  cross  -section  for  collisions 
of  ith-type  particles  with  jth-type  particles. 
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Because  of  their  extremely  small  mass,  the  electrons 
make  essentially  no  contribution  to  the  viscosity  of  the 
plasma  as  a  whole.   Thus,  equation  (29)  may  be  expressed  as 

-1 

where  ^N^[Tf^(v^Qu-fJ    •*■  ViQi^j]    ,  (33) 


-i 
and  \=[/^(rJHQi.H    -^^i^c'L)]    .  (34) 


For  plasmas  where  the  electron  temperature  is 
considerably  greater  than  the  heavy  particle  temperature, 
Fay  (2)    has  further  suggested  that  the  thermal  conductivities 
of  the  heavier  and  lighter  particles  be  considered  sepa- 
rately.  His  final  results  are  as  follows: 


4-  -m^ 


(35) 


-i 


where  the  subscript  h  refers  to  the  heavy  ion  and  neutral 

particle  components,  and  the  subscript  e  refers  to  the 

lighter  electron  component.   C,,   denotes  the  specific  heat 

Pn  ^ 

of  the  neutral  (or  ion)  particles.   Q^  is  the  Coulomb  cross 
section  for  charged  particles  and  is  given  in  dimensional 
form  as 
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Q,  =  i.soo-'j  r'^^^-^7.  (37) 

Simplified  relations  for  determining  the  appropriate  momentum 
transfer  cross  sections  from  cross  section  information 
commonly  available  are  given  by  Graf  (16) . 

A  more  practical  set  of  formulae,  as  far  as  ease  of 
numerical  calculations  is  concerned,  is  given  by  Powers 
(40) .   He  gives  the  viscosity  and  thermal  conductivity 
relations  in  dimensional  form  for  pure  argon  as  follows: 


jU^j^0.3iUio)%  -S.0O2(io'')Tj^  ■i-i-073(/d-'^)%  ,    (38) 


and     i{^^0.1d(i>7jLL,j.  (39) 

For  determining  electron  and  argon  ion  viscosities 
and  thermal  conductivities,  the  following  equations  were 
given: 

:=  0.589(10'^)  5  y0n[l.Z4(w)  %  ^7?^   7*  (41) 


'i^.rrf^'^/^     r.         .    -i-  J'^'^/     ^y^ 


/U.  =  117(10  )  n/3rr[l.2^(J0j  ^  /?7e  V,  C^2 ) 
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4  =0.Zi8(ld'*)  r!y£n[l24.(l0^)r-%'^''].  (43) 


As  the  gas  of  a  free  jet  flow  leaves  the  discharge 
region,  particle  density  and  temperature  gradients  are 
established  within  the  jet.   Because  of  these  gradients 
particle  diffusion  proceeds  from  regions  of  higher  toward 
regions  of  lower  concentration  and  temperature.   The 
electrons,  because  of  their  comparatively  low  mass  and 
consequent  high  mobility,  tend  to  diffuse  faster  than  the 
ions  or  neutral  particles.   The  ions,  on  the  other  hand, 
tend  to  remain  behind,  thus  creating  a  space  charge 
region  and  an  electric  field.   A  coupled  diffusion  process, 
known  as  ambipolar  diffusion,  results  in  which  the  electron 
diffusion  is  decelerated  and  ion  diffusion  is  accelerated. 
Both  particle  types  eventually  diffuse  with  the  same 
velocity . 

The  particle  density  level  above  which  space  charge 
effects  produced  by  the  interaction  between  electrons  and 
positive  ions  becomes  important  rather  than  independent 
electron  and  ion  diffusion  is  given  by  the  criterion  (41) 

'/&  >^ ~T^ — ~ •  (44) 

£q  is  the  permittivity  of  free  space,  e  is  the  electronic 
charge,  and  y  is  a  length  which  is  characteristic  of  the 
diffusion  region.   For  the  case  of  free  jet  mixing,  an 
appropriate  value  of  y  is  the  boundary  layer  thickness 
across  which  diffusion  takes  place. 

The  equations  governing  the  diffusion  of  electrons 
and  ions  in  the  neutral  parent  gas  of  a  weakly  ionized  flow 
are  given  by  (33)  and  (35)  as 
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fi  =  '^i^  =  -2)c^^i-*-y"i^s^i>  (45) 


and       #e  =  r7£  ^1  =  -D^  V??^  +/-<e  C  ^e  >  <  ^6 ) 


where  ^=  nv  is  the  mean  particle  flux,  e"^    is  the  electric 
field  created  by  the  space  charge  effects,  D  is  the  self- 
diffusion  coefficient,  and  u.     is  the  mobility  of  a  charged 
particle  as  defined  by 


V 

(47) 


The  subscripts  e  and  i  refer,  of  course,  to  the  electrons 
and  ions,  respectively. 

Equations  (45)  and  (46)  may  then  be  rewritten  as 


^  ="-^^^'^i-^A^5»  (48) 


and      %  ==-  ~Vy?e-t-JU^c^.  (49) 


If  the  charge  density  in  the  plasma  is  sufficiently 
high  for  ambipolar  diffusion  to  occur,  a  quasi-neutrality 
condition  ng  —  ni  then  exists.   Thus,  the  net  current  is 
zero,  and  n^v^  C;i  "e"^e  •   Thus, 


V-  -  Ve  =  - ..  (50) 
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The  coefficient 


is  known  as  the  ambipolar  diffusion  coefficient 


2.6   Langmuir  Probe  Theory 

One  of  the  earliest  and  most  useful  methods  employing 
probe  measurements  for  plasma  diagnostics  was  that 
employed  by  Langmuir  and  Mott-Smith  in  1924  (42) .   The 
Langmuir  single  probe  concept  has  been  found  useful  in 
obtaining  local  measurements  of  electron  temperature, 
electron  number  density,  and  wall  and  space  potentials. 
Although  the  "classical"  Langmuir  probe  theory  is  restricted 
to  stationary  and  particular  types  of  time-varying  discharges, 
it  is  probably  the  plasma  research  tool  which  has  provided 
the  bulk  of  experimental  knowledge  to  date.   Many  investi- 
gators have  used  this  single  probe  technique  and  found  the 
quantities  measured  to  compare  quite  favorably  with  those 
obtained  by  other  means  (1_6)  (43)  (£4)  (45)  .   Slight 
modifications  for  using  the  probe  techniques  in  circum- 
stances which  are  outside  the  limits  of  the  restrictions 
of  classical  theory,  e.g.,  in  a  strong  magnetic  field,  have 
also  been  developed  (46)  (47)  (48) . 

The  basic  theory  of  Langmuir  probes  has  been  outlined 
in  many  writings  on  the  subject.   Extensive  treatments  of 
this  theory  are  presented  in  references  (48)  (49)  (50)  (51). 
Therefore,  only  a  brief  summary  of  the  Langmuir  probe  and 
its  operation  as  pertinent  to  this  investigation  will  be 
included  herein. 
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In  the  Langmuir  single  probe  method  a  metallic  probe 
is  immersed  in  a  plasma.   The  probe  is  insulated  except  for 
a  small  area  which  is  exposed  to  the  plasma.   An  electric 
circuit  designed  so  that  the  probe  potential  may  be  varied 
with  respect  to  a  reference  electrode,  permits  the  measure- 
ment of  probe  current  as  a  function  of  its  potential. 
Depending  upon  the  voltage  impressed,  the  probe  will  either 
emit  or  collect  current.   An  analytical  interpretation  of 
the  current  drawn  by  the  probe  as  a  function  of  the  voltage 
applied  to  it  provides  information  about  the  local  conditions 
in  the  plasma.   The  current-voltage  curve  is  usually 
referred  to  as  the  probe  characteristic. 

The  basic  assumptions  of  classical  probe  theory  are 
outlined  by  Schwartz  (48)  as  follows: 

1.  The  sheath  thickness  is  small  compared  to  any  of 
the  mean  free  paths  in  the  system  under  considera- 
tion.  This  permits  the  region  of  effective  charge 
collection  to  be  considered  collisionless . 

2.  The  probe  diameter  is  small  compared  to  the  mean 
free  paths  in  the  system. 

3.  The  sheath  thickness  must  be  less  than  the  probe 
diameter . 

4.  A  Maxwellian  distribution  of  electrons  exists  at 
the  sheath  edge. 

If  the  electron  motion  in  the  plasma  is  purely  random, 
as  would  be  the  case  under  the  last  assumption  above,  the 
electron  concentration  at  the  probe  surface  n^^  is  given  by 
the  kinetic  relation 


^es=  i^ep  ^^P-ff^j*  (52) 
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where  n^   is  the  electron  concentration  outside  the  plasma 

sheath.   As  pointed  out  by  Gaither  (1^),  when  small  uniform 

field  forces  are  present,  this  equation  differs  only  by  an 

additive  constant. 

The  random  electron  current  density  to  the  probe  i„ 

«-s 

may  then  be  expressed  as 


k^K^p-Cw 


(53) 


where  the  random  plasma  current  density  j_   is  given  by  the 

*^P 

kinetic  relation 


K=    4.  (54) 


By  adding  the  ion  current  density  to  the  probe  i-:   to 

s 

the  electron  current  density,  the  total  random  current 
density  to  the  probe  can  be  obtained.   However,  in  the 
region  of  the  probe  characteristic  which  is  of  interest 
(the  retarding  field  region)  the  ion  current  drawn  by  the 
probe  is  negligible  compared  to  the  electron  current.   Thus 
ji   will  be  neglected  in  further  discussion. 

Taking  the  natural  logarithm  of  each  side  of  equation 
(53),  one  obtains 


Noting  that  ln(je  )  is  a  constant  and  differentiating  both 

Jr 

sides  of  this  equation  with  respect  to  ¥„  yields 


A(ill±^   =  _  _^_.  (56) 


Thus,  if  a  graphical  plot  is  made  of  ln(jg  )  versus  Vp  for 
the  retarding  field  region,  the  resulting  curve  should  be 
a  straight  line  with  slope  equal  to  -e/kT^.   If  the  curve 
is  indeed  a  straight  line,  the  original  assumption  of  a 
Maxwellian  distribution  is  verified  according  to  theory. 
If  the  curve  is  not  a  straight  line,  the  distribution  is 
either  non-Maxwellian  or  other  factors  have  affected  the 
resulting  curve,  e.g.,  an  unclean  probe  (49)  (50)  (51). 
The  electron  temperature  can  then  be  obtained  from 
the  equation 


Te=i.l<^(io)4S-7-^  (57) 


where  Tg  is  in  degrees  K  and  Vp  in  volts. 

From  kinetic  theory  the  average  electron  speed  in  a 
plasma  is  given  by 


The  electron  current  density  in  the  plasma  j^   may  also  be 
expressed  in  terms  of  the  electron  current  in  the  plasma 
io   and  the  probe  surface  area  A  as 
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By  combining  equations  (57),  (58),  and  (59),  the  following 
relation  is  obtained  for  calculating  the  electron  number 
density : 


^e^^^ie^i     (10  ),  (60) 


2 

where  A  is  in  cm  ,  ig   is  in  milliamps,  and  n^   is  in 

electrons/cm*^.   The  value  of  i^:,   to  be  used  in  this 

ep 

equation  is  that  corresponding  to  the  plasma  potential, 

i.e.,  at  the  break  in  the  curve  of  ln(i^  )  versus  V^  (19). 

tip  p  — 

Certain  criticisms  have  been  levied  at  the  use  of 
probes  to  measure  plasma  properties.   Chief  among  these  is 
that  the  insertion  of  the  probe  into  the  plasma  alters  the 
properties  to  be  measured.   Other  criticisms  pertain  to  the 
reflection  and  secondary  emission  of  electrons  by  the  probe 
and  to  possible  thermal  electron  emission.   Extensive 
discussions  of  these  possible  sources  of  error  are  given 
in  references  (48)  (49)  (50)  (51). 


3.   THEORY  OF  THE  EXPERIMENT 

This  study  is  a  part  of  a  large  scale  effort  to  obtain 
a  complete  solution  and  understanding  of  the  problem  of 
charged  particle  behavior  in  low  pressure  plasma  flows.   The 
complete  problem  involves  simultaneous  consideration  of  the 
macroscopic  character  of  the  flow  as  well  as  individual 
component  behavior.   Fortunately,  in  flow  systems  such  as 
the  one  under  consideration  for  cases  of  weak  ionization, 
the  macroscopic  properties  are  affected  very  little  by  the 
fact  that  the  gas  is  slightly  ionized.   Thus,  these  proper- 
ties may  be  examined  with  the  knowledge  that  the  macroscopic 
behavior  of  the  jet  flow  differs  very  little  from  that  of 
an  unheated  gas  flow. 

Earlier  studies  of  the  process  of  free  expansion  of  a 
gas  from  a  nozzle  or  orifice  have  been  made  by  several 
investigators  using  both  macroscopic  and  microscopic  view- 
points.  Pai  (2^)  (11)  studied  the  laminar  mixing  of  neutral 
gas  flows  and  was  able  to  obtain  complete  analytical  solu- 
tions for  the  macroscopic  velocity  and  temperature  profiles. 
Gaither  (15)  and  Graf  (16)  investigated  free  expansion 
argon  plasma  jets  and  were  able  to  determine  the  charged 
particle  behavior  in  the  inviscid  core  and  developed 
regions,  respectively.   Thus,  the  analytical  description 
of  flr,vv  processes  in  the  region  near  the  nozzle  is  essen- 
tially complete  except  for  knowledge  of  charged  particle 
behavior  in  the  jet  mixing  and  base  regions.   It  is  this 
author's  intent  to  complete  this  picture. 
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In  this  study  considerable  effort  was  expended  in  the 
control  and  elimination  of  complexities  that  might  be  found 
in  flow  systems  outside  the  laboratory.   However,  the 
complex  character  of  the  flow  was  still  such  that  certain 
simplifications  had  to  be  made  to  obtain  meaningful  infor- 
mation.  The  objective  of  the  work  was  accomplished  using 
a  relatively  simple  theoretical  model  having  the  following 
restrictions: 

1.  Only  three  species  need  to  be  considered:  (a)  argon 
atoms  in  the  ground  state,  (b)  singly  ionized  argon 
ions  in  the  ground  state,  and  (c)  electrons. 

2.  Each  component,  as  well  as  the  mixture,  obeys  the 
perfect  gas  equation  of  state. 

3.  Local  thermodynamic  equilibrium  is  maintained 
throughout  the  jet  flow. 

4.  Radiation  effects  are  negligible. 

5.  The  flow  is  axisymmetric . 

6.  The  flow  is  steady. 

7.  The  flow  is  laminar. 

8.  Compressibility  effects  are  negligible. 

9.  The  pressure  is  constant  throughout  the  region  of 
interest . 

10.  Specific  heats  are  constant. 

11.  The  boundary  layer  approximation  is  valid  for  the 
mixing  zone. 

12.  The  average  properties  of  individual  components 
are  appropriately  described  by  classical  mechanics 
(see  Appendix  B) . 

13.  Charged  particle  diffusion  takes  place  as  ambipolar 
diffusion . 

14.  There  are  no  externally  applied  fields  (including 
gravity),  except  that  of  the  base  region  which  is 
considered  separately. 


The  validity  for  establishing  these  restrictions  has  either 
been  brought  out  or  will  be  discussed  later  in  connection 
with  the  experimental  results . 

In  keeping  with  the  aforementioned  boundary  layer 
approximation  as  it  applied  to  a  circular  jet  mixing  region, 
the  equations  of  continuity,  momentum,  and  energy  for  the 
axially  symmetric  laminar  mixing  of  an  incompressible  fluid 
may  be  written: 


u  and  V  are  the  velocity  components  in  the  z  and  r  directions, 
respectively.    u.  is  the  viscosity,  and  K"  is  the  thermal 
conductivity. 

With  no  pressure  gradients  in  the  flow,  a  Prandtl 
number  equal  to  unity,  and  constant  specific  heats,  Pai  (11) 
shows  that 


T=A-»-6u.-  a/aCp. 


(64) 


T  is  the  temperature  at  a  given  point  in  the  jet,  and  u  is 
the  corresponding  axial  velocity.   A  and  B  are  constants  to 
be  determined  by  the  boundary  conditions. 

For  the  case  where  the  velocity  and  temperature  in  the 
jet  differ  only  slightly  from  the  surrounding  quiescent 
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medium,  Pai  (11)  uses  small  perturbation  techniques  and 
the  boundary  conditions  at  the  nozzle  exit,  namely 

at  z  =  0:    u  =  u^  =  constant  for  0  —  r  ~  Rq 

u  =  Uq  for  r  >  Rq 

to  obtain  the  velocity  distribution  relation 


^=\/  e  JM)J,m:)<^'^'  (65) 


The  subscript  0  refers  to  some  reference  state,  and  7\     is 
the  eigenvalue  determined  by  the  boundary  conditions.   R 
is  the  nozzle  exit  radius,   Jq  and  S-^   are  the  bessel  func- 
tions of  zeroth  and  first  order  of  the  first  kind,  respec- 
tively.  Thus,  the  velocity  distribution  may  be  plotted 
for  various  values  of   c^'^z,  where   cx   is  defined  by 


^   Mo 

^'^-KTT'  (66) 


PoUo 


By  introducing  the  stream  function  into  the  equation  of 
motion,  Pai  (1_1^  takes  a  more  rigorous  approach  to  obtain 
an  exact  velocity  distribution  solution  for  the  laminar 
case.   His  final  equation  can  be  integrated  numerically 
using  a  step-wise  finite  difference  method  to  give  the 
velocity  distribution.   The  results  of  this  method  and  that 
of  the  small  perturbation  method  are  shown  in  graphical 
form  for  various  values  of  o<  ^x  in  reference  (11), 

Of  fundamental  concern  in  this  investigation  are  the 
spatial  variations  of  electron  number  density  and  electron 
temperature  in  the  jet  mixing  and  base  regions  of  such  a 
jet  flow  field.   An  extensive  investigation  of  these 
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parameters  in  the  inviscid  core  region  of  a  low  pressure 
argon  plasma  jet  has  been  made  by  Gaither  (15) .   He  found 
the  electron  temperature  in  this  region  to  be  constant  in 
the  radial  direction.   In  the  axial  direction  the  electron 
temperature  dropped  off  in  an  exponential  manner.   The 
electron  number  density  in  the  radial  direction  correlated 
quite  well  with  his  theoretical  prediction  of  behavior  as  a 
Bessel  function  of  zero  order.   It  also  dropped  off  expo- 
nentially in  the  axial  direction. 

Graf  (16)  investigated  a  low  density,  free-expansion 
argon  plasma  jet  and  examined  the  electron  density  profiles 
in  the  developed  region  downstream  of  the  core  region.   The 
electron  density  was  found  to  have  a  Gaussian  radial  distri- 
bution, while  falling  off  approximately  exponentially  along 
the  jet  axis.   No  information  regarding  the  electron 
temperature  profiles  in  the  jet  was  given.   Neither  Graf 
nor  Gaither  relied  upon  Pai's  theoretical  model  development. 

Radial  electron  temperature  profiles  in  the  positive 
columns  of  low  pressure  discharge  tubes  have  also  been  found 
to  be  essentially  constant  (33)  (35) .   The  question  then 
arises  as  to  whether  or  not  such  a  characteristic  is  also 
maintained  throughout  the  base  and  boundary  layer  regions 
in  a  jet  flow  field.   In  particular,  it  is  of  interest  to 
know  if  the  electron  temperature  also  remains  constant 
radially  throughout  the  region  in  which  the  gas  mixes  with 
the  surrounding  medium.   Perhaps  the  mixing  action  of  this 
region  affects  the  electron  behavior  in  such  a  manner  as  to 
cause  the  electron  temperature  to  relax  toward  the  tempera- 
ture of  the  parent  gas. 

Cobine  (19)  has  shown  that  for  a  given  gas  the  following- 
relation  holds  for  the  positive  column: 


v.-  Vji  1  1  itf  ■> 


whore  P  is  the  pressu:.'.      .^    \...^   -..^^^^._.   ^^  .. 
xube,  dAC   C  is  u  constai-^  v.aich  depends  on  the  ^as  (4  x  10"*=^ 
for  argon)  ,  CJ  is  defir.     .  terms  of  the  ioni  .ation 
potential  Vi  as 

Thus,  the  electron  temper:  .iction  only  of  pres- 

sure and  tube  radius  in  the  positive  colurr.n  of  a  .':':iven  gas 
dischar^^.   rhe  possibili-^  ^..-...  ^.„^.:^  „..^.  ...^  _.^...^ 
characteristic  will  be  displayed  by  the  flow  upo.". 
the  discharge  region. 

A  theoretical  evaluation  ^: 

» 

electron  density  distribution  : 

nozzle  exit  can  bo  r.ado  by  con^ 

vatioi:  equation  -v^_  \...^    -._w.  „„^^..  :._  c.-.^,. 

reference  (15)  in  the  lorn 


where  q  is  the  net  electron  velocity  due  to 
forces,  and  N  is  the  net  volume  rate  of  electi 
D^  and  Dr  are  the  ambipolc.i->  diffusior,  coeffic" 
due  'cw  ...^-  .v^rticle  der.o_cj  and  t_....   ature  ^- 
within  -che  jet  flov/.   Equation  i^'^^)    can  be  : 
form 


(70) 


if  the  1 '.-llov/in':;"  criteria  are  met: 


d  recombinatic.  o^x\;^„o 


icantly  sinall,  i.e.,  K  ^  0, 

■q  may  be  adequately  descri". 

velocity  q^ . 

The   ai.. bipolar   diffusion   coc  i:'icients   are 


4.  The  radiul  diffusion  current  Cv;h.     ^curs  duo  to 
the  cc,.»,^^..w>..  effects  of  spev-.^>-.o  ^j.c<v-j.^.ii  .=,  c^.,v. 
species  energy  gradient)  is  large  coiVxpared  to 
axial  diffusion  current. 

5.  Tg  is  essentially  constant  in  the  radial 

The  boundary  conditions  cw  \j^    ^.,.^^^^<..   w..  ov_....  .ion  (70) 

are  those  at  the  nozzle  mouth  anc  the  test  section  walls. 

Since  in  this  investigation  the  ,r   is 

approximately  10  times  the  noz:-!.  region 

of  interest  is  confi/.ed  to  tha":  ^t 

flow  may  be  a:_^. .„  „„  .-„„.....^  _..^^  „..  ...-.-.. 

space  region.  Tiie   boundary  conditions  to  be  i:r.pct^ww  :..;^;,' 

then  be  expressed  as 

i:  r.   =    finite  -    0: 

ng   =   0  J  r  ^  11^,,    z   =   0; 

r-^   =  ne(-y  ^--  r  •<.  Rq,    z   =   0;                                      (71) 

-^ — ^0  as  r — 'j^^-L.  ; 

iiQ — --0  as  z—i^-c.   . 

Ey  assuming  n^  =  ^  yj.  j ..-i  ^r)  ,    z.ia   me'cnoc  ^^^    separation 
of  variables  may  be  imple/r.ented,  and  the  solution  (52) 
zo   equation  (70)  may  be  expressed  as 


This  :ivjiu.;ion  i.:ay  oe  presentee  in  a  torm  wnicL  :^s  part 
larly  useful  in  this  study  as 


The  u.^o  oi  an  i..     -i  lorm  on  tn. 
of  a  sui.imatic     ..pplicable  if  the 
pass  from  the  nozzle  into  an  unboi. 
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spacing  bet\yeen  the  eigenvalues  decr^_-..^..  «..  ^.,^    ^^.^^    ^^ 
the  boundaries  increases  (_54)  . 
Define  f(r)  such  that 


T^'O  ^  '  /c-  i't-'y-'  ''• '  ''■'  *^-  i  •'• ''/ 


/;  -.  <• . 
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The  A(A)  values  may  be  cleterinined  by  multiplying  both  side 
of  equation  (74)  by  r  Jq(^  r)dr  and  integrating  from  0  tc 


rvO  i-j  -^  -  'V  ■'- ■  "        -  y  .  I  - .-  .       ( 75 ) 


From  xhe  Fourier-Bessel  integr 


Therefore , 


since  ^   can  be  cho^     .ual  to       Thus 
function  f  (r)  tho  values  of  A(^\  )  can  be  :_ 

tion  of  A(?v)  ^...w  ^v-^c...  o-wn  (73;.   -^_v..o  ...^  

for  the  number  density. 

The  fuiiction  f(r)  is,  of  course- 
density  profile  at  the  nozzle  exit.   This  .for.T.a- 
tion  may  be  obtained  experimentally,  or  a 
assumed,   irov/cvc:-.  the  closer  to  "Jr..    r.oz~" 
density  is      -  determined,  thv^  ...w.^  ...^^ .  .  ..^  ...^  ;^.,:i. ,_ 
for  f(r: 


The  initial  electron  number  density  profij.^ 
plasma  jet  has  oee-A   shown  to  be  appropi-iately  descrx^ea, 
both  theoretically  and  experimentally,  by  a  Bessel  function 
of  zero  order  and  first  kind  (15)  (19).   However,  all  pre- 
vious theoretical  analyses  have  been  based  on  the  assumption 
that  yiq   approached  zero  at  r  =  0  for  all  values  of  z . 
assumption  is  not  entirely  cori-^^o,  ^:^..c^-  :.  c  does  not  .._:..^.. 
for  radial  spread  of  the  electron  nu.aber  density  profiles 
as  the  jet  proceeds  from  the  nozzle.   The  approach  taken 
herein  eliminates  that  assumption. 


4.   EXPERIMENTAL  CONSIDERATIONS 

4.1   Description  of  Experimental  Apparatus 

The  experiments  were  carried  out  in  the  Graduate  Research 
Laboratory  of  the  Mechanical  Engineering  Department  at  the 
University  of  Florida.   The  general  features  of  the  experi- 
mental apparatus  and  related  instrumentation  are  shown 
schematically  in  Figure  3.   The  test  gas  was  supplied  from 
a  high-pressure  bottle  through  a  flow  regulator  and  control 
valve  to  a  settling  chamber  before  entering  the  plasma 
generator  proper,  where  it  was  ionized.   Argon  was  chosen 
as  the  working  gas  since  it  was  free  from  any  tendency  to 
form  negative  ions  upon  being  ionized  and  its  physical 
properties  were  well  known.   The  ionized  gas,  or  plasma, 
then  expanded  through  a  converging  nozzle  into  the  test 
section.   A  continuously  operating  mechanical  vacuum  pump 
then  removed  the  gas  from  the  test  section.   Figures  4  and  5 
show  the  facility  and  test  apparatus  used  for  the  experiment. 

The  following  is  a  description  of  the  apparatus  wherein 
each  of  the  components  is  explained  separately. 

Plasma  Generator 

The  plasma  generator,  shown  in  Figure  5  and  schematically 
in  Figure  1,  consisted  of  five  basic  components:  the  settling 
chamber  consisting  of  several  baffling  chambers,  the  chamber 
nozzle,  the  discharge  zone  housing,  the  cathode  rod,  and  the 
anode  plate  assembly.  The  gas  initially  entered  the  base  of 
the  generator  through  two  1/4  inch  inlets  and  then  proceeded 
to  flow  through  the  baffling  chambers,  which  served  to  direct 
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the  flow  uniformly  and  coa.xially  through  the  chamber  nozzle. 
This  nozzle,  constructed  of  lucite,  further  helped  to 
the  gas  flow  unifoi-m  upon  entering  the  dischar         i . 

A  1/4  inch  OD  thin-walled  stainless  steel  tu„_ 

as  -Che  cathode.   One  end  of  the  tube  was  v/elded  and  grounded 
to  a  sr.iooth  parabolic  cor.-cour.   Through  tl.e  open  end  of  the 
cathode  a  1/3  inch  OD  tube  v/as  passed  which  allov/ed  the 
cathode  to  be  air  cooled.   The  cathode  was  designed  such 
that  it  was  vacuum  sealed  using  o-rings  and  could  be 
positioned  ax  any  desired  distance  ircr.  the  ano^c 

The  anode  plate  assembly  c.,...^_^  .^.,  ^_  . 

i/2  inch  thick,  ,_ ,-rt  plate  --co   whic" 

inserted.   The  co^iVerging  nozzle  was  cons - 

less  steel  and  was  designed  so  zL:.z    it  could  be  c. 

necessary.   Experiments  showed  that  a  nozzle  exit  d; 

of  15/32  inch  (11.91  mm)  was  suitable  for  the  inv- 

Reasons  for  choosing  this  partic.._._  size  nozzle  are  dis- 
cussed later  in  this  chapter. 

The  discharge  zone  housing  was  construe        .  75  mm 
OD  Pyrox  glass  tube,  15  cm  in  long...   ....  :...  •  ;,,  .  .^  qd 

Pyrex  glass  tube  was  added  after  the      .al  f^enerator 

construction  to  help  cli-:inate  . .._ _..  ^._  . ^v...^.-^\ 

zone.   The  transparent  .  _  .-ex  glass  housing  provide    "  .- c- 
trical  insulation  and  permitted  better  adjustment  of  che 
gas  flow  rate  to  yield  a  plasma  jet  that  was  as  1..        ; 
possible.   Vacuum  sealing  was  obtained  by  using  1/32  inch 
thick  neoprene  gaskets  on  the  ends  of  the    .  -j . 

Base  Region  At...  ....  , 


The  base  ro  .Ion  assembly  was  actual!,,- 
Side  of  the  z...  .  te   assembly  she, 

to  the  i.acure  of  the  investigation 
able  to  va.-      potential  in  this  . 
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the  anode  and/or  the  exhaust  lead  housing  to  see  its  effect 
upon  the  electron  diffusion  process.   Thus,  it  was  necessary 
to  isolate  this  region  from  the  anode  plate  assembly.   This 
isolation  was  obtained  by  placing  a  1/16  inch  thick  lucite 
disk  between  a  1/16  inch  thick  copper  disk  and  the  anode 
assembly . 

Test  Section  and  Vacuum  System 

The  test  section  was  a  vertically  mounted  100  mm  ID 
Pyrex  glass  tube,  75  cm  long,  supported  by  and  sealed  to 
the  anode  assembly  and  the  exhaust  lead  housing.   Static 
pressure  taps  were  located  along  the  side  of  the  test 
section  at  distances  of  3,7  cm,  7.8  cm,  11.7  cm,  and  19.7 
cm  from  the  base  region. 

A  Consolidated  Vacuum  mechanical  vacuum  pump,  with  a 
225  m^/hr  pumping  capacity,  maintained  a  pressure  of  40-60 
microns  Hg  in  the  test  section  when  the  test  gas  flow  was 
shut  off.   A  nominal  mass  flow  rate  of  0.2  gf^/sec  raised  the 
pressure  to  3-4  mm  Hg.   A  surge  control  valve  located  in 
the  exhaust  line  to  the  vacuum  pump  and  a  valve  in  the 
exhaust  lead  housing,  permitting  test  gas  to  be  introduced 
downstream  of  the  plasma  jet,  allowed  the  back  pressure  in 
the  test  section  to  be  adjusted  quite  accurately. 

Electrical  System 

Two  power  supplies  of  the  type  needed  to  create  the 
discharge  were  available  for  the  experiments:  a  0-3  ampere, 
0-1500  volt  DC  Rapid  Electric  Silicon  Rectifier  and  a  0-1 
ampere,  0-2000  volt  DC  High  Voltage  Generator  Unit.   Both 
power  supplies  were  used  in  the  preliminary  testing  and 
neither  showed  any  particular  advantage  over  the  other.   The 
latter  was  chosen  for  use  with  the  plasma  generator  in  the 
final  experiments.   To  assist  in  current  and  voltage  control 
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of  the  electrical  discharge,  a  0-2040  ohm  High  Voltage 
Ballast  Resistor  Unit  was  connected  in  series  to  the  cathode 
side  of  the  power  supply.   The  anode  plate  assembly  was 
grounded.   Pi'obe  circuitry  is  discussed  later  in  this  section. 

Cooling  System 

Initial  system  design  provided  for  internal  cooling  of 
both  the  anode  and  cathode  assemblies, using  laboratory 
supplied  compressed  air.   Electrical  insulation  and  flexi- 
bility were  obtained  by  using  1/4  inch  high  pressure  rubber 
hoses,  fitted  with  1/4  inch  rubber  hose  fittings.   Early 
experiments  indicated  that  it  was  unnecessary  to  cool  the 
anode  plate  assembly  due  to  the  low  degree  of  ionization  in 
the  plasma  jet.   A  rather  large  volume  of  air  was  required, 
however,  to  properly  cool  the  cathode. 

Probe  Positioning  Mechanism 

Design  of  the  exhaust  lead  housing  provided  for  the 
insertion  of  three  probes  into  the  flow  stream.   These 
probes  were  positioned  by  the  probe  positioning  mechanism. 
This  device  was  located  on  top  of,  but  electrically  insulated 
from,  the  exhaust  lead  housing. 

All  positioning  of  the  probes  was  accomplished  manually. 
Gear  ratios  were  such  that  the  probes  could  be  positioned 
axially  to  the  nearest  0.30  mm  and  radially  to  the  nearest 
0 . 06  mm. 

Probes  and  Instrumentation 

The  plasma  jet  flow  field,  shown  in  Figure  6,  was 
surveyed  using  three  types  of  probes:  an  impact  pressure 
probe,  a  thermocouple  probe,  and  a  Langmuir  probe.   All  three 
probe  stems  were  constructed  of  5/16  inch  OD  Mulltite 
ceramic  tubing  which  provided  the  necessary  electrical 
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insulation.   A  1  3/8  inch  brass  offset  was  attached  to 
each  probe  stem,  permitting  the  probe  to  be  moved  radially- 
through  the  plasma  jet  as  the  stem  was  turned.   Wilson -type 
seals  were  used  in  the  exhaust  lead  housing  to  insure  proper 
vacuum  sealing  of  the  probes. 

Impact  pressures  were  measured  with  an  open-ended 
brass  tube  (attached  to  the  brass  offset)  which  pointed 
into  the  flow  stream.   The  tube  had  an  inside  diameter  of 
0.067  inches  and  was  cut  to  a  20  degree  chamfer.   The  other 
end  of  the  probe  was  connected  with  Tygon  vacuum  hose  to  a 
0.1-20  mm  Hg  Wallace  Si  Tiernan  Absolute  Pressure  Indicator. 
Corresponding  static  pressures  were  obtained  by  recording 
the  pressure  when  the  impact  pressure  probe  was  completely 
outside  the  jet  flow. 

Stagnation  temperatures  were  measured  using  a  shielded 
thermocouple  constructed  of  0.010  inch  diameter  chromel  and 
alumel  wires.   A  thermo  Electric  Potentiometer  Pyrometer 
was  used  to  obtain  the  temperature  readings  directly  in 
degrees  Farenheit. 

Both  cylindrical  and  flat  Langmuir  probes  were  used  for 
the  investigation.   The  plane  probes  were  constructed  of 
0.020  and  0.030  inch  diameter  tungsten  wire,  while  the 
cylindrical  probes  were  constructed  of  0.005  inch  diameter 
wire  with  the  exposed  surface  being  approximately  0.065  inches 
in  length.   The  Pyrex  glass  insulation  on  the  probe  was  made 
as  small  as  possible  so  as  not  to  disturb  the  flow  stream. 
Since  it  was  difficult  to  clean  the  collecting  surface  near 
the  insulation,  an  attempt  was  made  to  keep  the  glass  from 
touching  the  probe  metal  in  the  vicinity  of  the  collecting 
surface.   Figure  7  shows  a  sketch  of  the  Langmuir  probe 
construction. 

The  voltage-current  characteristics  for  the  Langmuir 
probe  were  obtained  using  the  circuit  shown  in  Figure  S.   A 
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110  volt,  60  cps  Variac  was  used  to  establish  the  desired 
potential  at  point  A,   A  one-to-one  isolation  transformer 
was  used  purely  for  isolation  purposes.   By  suitable  adjust- 
ment of  the  Variac,  resistor  R-l,  and  the  probe  bias,  the 
probe  potential  Hb  could  be  swept  over  any  desired  amplitude 
about  a  suitable  bias.   The  probe  potential  was  input  to  the 
horizontal  deflection  plates  of  a  Tektronix  Type  502  Dual- 
Beam  Oscilloscope. 

The  probe  current  was  measured  across  resistor  R, , 
Both  resistors  R-^   and  R2  were  made  variable  so  as  to  give 
flexibility  in  measuring  the  electron  current  at  positions 
throughout  the  plasma  jet,  although  R^  had  to  be  accurately 
measured  using  a  resistance  box  for  each  probe  reading. 
Since  the  signal  V^;-  -  Vg  was  proportional  to  the  current 
drawn  across  R-j^,  the  voltages  V^-;  and  Vg  were  input  to  a 
differential  amplifier  on  the  vertical  deflection  plates  of 
the  oscilloscope. 

The  maximum  current  across  R-,  was  approximately  10 
milliamps  near  the  centerline  of  the  jet.   Near  the  edge  of 
the  jet  the  current  was  considerably  less.   Thus,  resistor 
R2  was  necessary  to  account  for  the  finite  input  impedance 
(1  Meg  ohm)  of  the  oscilloscope.   The  proper  setting  for  R2 
was  obtained  by  disconnecting  the  probe  from  the  probe 
circuit  at  the  BNC  connector  on  the  circuit  box  and  observing 
the  oscilloscope  trace.   It  was  desired  to  have  a  horizontal 
line  as  the  trace  under  such  conditions,  and  this  was 
usually  the  case.   However,  for  large  values  of  R-^   and  high 
vertical  displacement  sensitivities  the  trace  was  sometimes 
inclined.   The  resistor  R2  was  then  adjusted  so  as  to  make 
the  trace  become  horizontal  and  thus  form  a  consistent  base 
line  from  which  to  measure  the  vertical  deflections. 

A  typical  Langmuir  probe  characteristic  is  shown  in 
Figure  9.   The  corresponding  current-voltage  curve  is  plotted 
in  Figure  10. 
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4.2   Preliminary  Testing 

A  rather  extensive  preliminary  investigation  was  made 
to  determine  which  plasma  generator  configuration  would 
yield  plasma  flows  suitable  for  studying  jet  mixing  and 
diffusion  processes.   This  required  the  production  of  a 
stable  jet  large  enough  in  size  to  permit  insertion  of 
conventional  probes  and  yet  sufficiently  ionized  that  the 
plasma  properties  could  be  determined. 

Initial  generator  designs  allowed  for  air  cooling  of 
both  cathode  and  anode  assemblies.   Functional  tests  showed 
that  it  was  not  necessary  to  cool  the  anode  assembly.   Thus, 
this  assembly  was  reduced  in  relative  physical  size  to  that 
shown  in  Figure  1.   Using  air  to  cool  the  cathode  presented 
two  rather  complicated  problems:  (a)  a  large  amount  of  drift 
in  the  electron  temperature  and  electron  density  in  the 
plasma  jet,  and  (b)  the  appearance  of  dark  blue  spots  on 
the  cathode  from  time  to  time,  causing  the  Langmuir  probe 
characteristics  to  become  very  erratic.   The  first  problem 
was  attributed  to  the  fact  that  the  air  compressor  used 
was  positioned  too  close  to  the  test  apparatus.   This 
permitted  all  pressure  variations  in  the  air  storage  tank 
to  be  sensed  by  the  cathode.   The  blue  spots  which  occasion- 
ally appeared  on  the  cathode  were  attributed  to  a  slight 
condensation  of  the  water  vapor  in  the  cooling  air, 
resulting  in  an  uneven  cooling  of  the  cathode.   A  filter 
trap  was  installed  in  the  air  line,  essentially  eliminating 
these  two  problems.   An  isolated  water  supply,  employing  a 
recirculating  pump,  was  also  built  and  found  to  be  rather 
effective  in  eliminating  these  cooling  problems.   However, 
it  appeared  that  either  a  larger  cathode  or  a  higher 
pumping  pressure  would  also  be  needed. to  permit  water  cooling, 

It  was  desired  to  have  the  plasma  generator  operate 
over  as  wide  a  flow  and  pressure  range  as  possible.   Early 
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generator  designs  showed  very  turbulent  flow  conditions 
inside  the  discharge  housing  near  the  anode  plate  assembly. 
The  installation  of  a  small  Pyrex  glass  tube  inside  the 
discharge  housing  and  concentric  to  the  cathode,  as  shown 
in  Figure  1,  served  to  better  guide  the  flow  and  eliminate 
this  turbulence. 

Throughout  the  preliminary  investigation  numerous 
Langmuir  probe  surveys  of  the  plasma  jet  flow  field  were 
made.   Results  from  the  cylindrical  probes  compared  favor- 
ably with  the  plane  probe  results;  however,  the  break  in 
the  electron  saturation  part  of  the  probe  characteristics 
was  more  defined  for  the  cylindrical  probes.   For  this 
reason  and  due  to  the  fact  that  the  cylindrical  probes 
appeared  to  disturb  the  flow  less  than  the  plane  probes, 
they  were  chosen  for  use  in  the  final  testing. 

Converging  nozzles  of  1/4,  15/32,  and  9/16  inch  diam- 
eters were  used  to  determine  which  size  was  best  suited  for 
the  experiment.   The  1/4  inch  nozzle  did  not  permit  easy 
access  for  probing  the  flow  field  and  was  eliminated  for 
that  reason.   However,  it  did  yield  nozzle  exit  electron 
temperatures  on  the  order  of  8,000  to  12,000  X  which  were 
comparable  to  those  obtained  by  Gaither  (15)  under  quite 
similar  conditions.   Both  the  15/32  and  9/16  inch  diameter 
nozzles  provided  more  room  for  probe  measurements.   The 
15/32  inch  diameter  nozzle  was  chosen  since  under  the  same 
flow  conditions  it  would  yield  considerably  higher  electron 
temperatures. 

Total  pressure  surveys  for  each  of  the  above  nozzles 
indicated  the  presence  of  a  rather  thick  boundary  layer 
build-up  at  the  nozzle  exit.  Slight  changes  in  nozzle 
contours  had  negligible  effect  on  the  presence  of  the 
boundary  layers.   These  boundary  layers  will  bo  discussed 
in  more  detail  in  conjunction  with  the  final  test  results. 
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4.3  Experimental  Procedure 

Upon  completion  of  the  preliminary  testing,  a  series 
of  test  runs  was  made  to  determine  the  effects  of  pressure, 
power  input,  electrode  separation,  and  mass  flow  rate  on 
the  electron  temperature  and  electron  number  density 
distributions  in  the  jet.   Each  of  these  runs  was  given  a 
corresponding  test  number,  e.g.,  test  I,  test  II,  etc.   A 
complete  run  consisted  of  the  gathering  of  Langmuir  probe, 
pressure,  and  temperature  data  from  the  entire  flow  field 
up  through  the  initial  developed  flow  region. 

Each  run,  usually  lasting  from  six  to  eight  hours,  was 
begun  by  allowing  the  plasma  jet  to  stabilize  for  about  one 
hour.   After  the  stabilization  period  all  pertinent  pre-run 
information  was  recorded,  followed  by  the  recording  of  the 
Langmuir  probe  data.   These  probe  characteristics  were 
recorded  by  taking  oscilloscope  photographs  of  the  traces 
corresponding  to  the  various  probe  positons.   Experimenta- 
tion showed  that  several  traces  could  be  recorded  on  the 
same  photograph  by  superimposing  them.   Langmuir  probe 
readings  were  taken  at  radial  intervals  of  approximately 
0.3  mm  and  at  axial  intervals  of  approximately  10  mm.   Near 
the  edge  of  the  jet  the  radial  readings  were  sometimes  taken 
at  approximately  0.15  mm  intervals  in  an  attempt  to  better 
define  the  profile  information  in  the  mixing  region. 

Upon  completion  of  gathering  the  Langmuir  probe  data, 
impact  pressure  readings  were  taken,  followed  by  thermo- 
couple probe  readings.   Both  pressure  and  temperature 
readings  were  taken  at  radial  intervals  of  approximately 
0.6  mm  and  axial  intervals  of  approximately  10  mm.   The 
run  was  ended  by  the  recording  of  all  pertinent  post-run 
information. 

The  recording  of  the  necessary  Langmuir  probe  data 
usually  took  from  two  to  three  hours.   Slow  variations  in 
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electron  density  were  noted  during  this  time  by  making 
periodic  checks  near  the  jet  centerline  or  by  monitoring 
a  given  point  for  a  period  of  time.   These  variations  were 
usually  no  more  than  about  10%  of  the  original  reading  at 
the  point  under  consideration.   No  attempt  was  made  to  make 
any  adjustments  for  these  variations. 


5.   DISCUSSION  OF  EXPERIMENTAL  RESULTS  AND 
COMPARISON  WITH  THEORY 


Following  the  preliminary  testing  previously  outlined, 
several  test  runs  were  made  in  direct  accordance  with  the 
specific  objectives  of  the  investigation.   The  results  of 
three  of  the  tests  are  reported  herein.   The  controlled 
conditions  under  which  the  tests  were  conducted  are  as 
follows: 

Test  I 

Argon  Inlet  Temperature  300  K 

Discharge  Voltage  370  volts 

Discharge  Current  600  ma 

Jet  Static  Pressure  5.29  mm  Hg 

Mass  Flow  Rate  0.227  g^/sec 

Potential  Core  Velocity  220  m/sec 

Electrode  Separation  Distance  55.6  mm 

Base  Plate  Potential  (floating) 

Base  Plate  Current  0.0  ma 

Test  II 

Argon  Inlet  Temperature  298  K 

Discharge  Voltage  424  volts 

Discharge  Current  870  ma 

Jet  Static  Pressure  3.28  mm  Hg 

Mass  Flow  Rate  0.222  gn^/sec 

Potential  Core  Velocity  304  m/sec 

Electrode  Separation  Distance  55.6  mm 

Base  Plate  Potential            .  (floating) 

Base  Plate  Current  0 . 0  ma 
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Test  III 

The  test  conditions  for  this  run  were  the  same  as  those 
for  Test  II,  with  the  following  exceptions: 

Base  Plate  Potential  +71  volts 

Base  Plate  Current  145  ma 

The  particular  experimental  objectives  of  the  first 
two  tests  were  as  follows:   (a)  to  obtain  gas  dynamic  infor- 
mation in  order  to  define  flow  characteristics  and  identify 
different  flow  regions,  (b)  to  obtain  electron  temperature 
and  electron  number  density  profiles  throughout  the  entire 
flow  field  near  the  nozzle  and  as  far  away  radially  from 
the  jet  centerline  as  instrumentation  would  permit,  and 
(c)  to  examine  the  effects  of  certain  variables,  e.g., 
pressure  and  power  input,  on  these  profiles.   The  primary 
purpose  of  Test  III  was  to  determine  the  extent  to  which 
an  applied  potential  in  the  base  region  would  alter  the 
electron  diffusion  characteristics  of  the  jet  and  thereby 
establish  what  changes  in  electron  temperature  and  electron 
number  density  profiles  could  be  effected  by  this  means. 
The  basic  test  conditions  for  this  latter  run  were  inten- 
tionally kept  the  same  as  those  for  Test  II  for  purposes 
of  exact  comparison.   Experimental  results  for  these  tests 
are  presented  graphically  in  Figures  11  through  24. 

Test  I  was  run  at  an  exit  flow  Mach  number  (at  jet 
centerline)  of  0.62,  while  for  Test  II  the  Mach  number  was 
0.91.   Examination  of  the  corresponding  velocity  profiles 
in  Figures  11  and  16  reveal  the  presence  of  initial  boundary 
layers  at  the  nozzle  exit.   The  presence  of  such  thick 
boundary  layers  served  two  very  useful  purposes:   (a)  they 
gave  the  jet  mixing  region  an  initial  thickness,  thus 
making  it  wider  and  more  accessible  to  probe  studies,  and 
(b)  thoy  permitted  electron  density  and  electron  temperature 
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profiles  to  be  determined  within  this  initial  region  at 
the  nozzle  exit. 

Reasonable  justification  for  the  existence  of  these 
relatively  thick  boundary  layers  was  made  by  a  very 
simple  theoretical  analysis.   The  Reynolds  number  Rg, 
based  on  nozzle  diameter  and  exit  flow  conditions  at  the 
center line,  was  determined  to  be  980  and  1038  for  Teste  I 
and  II,  respectively.   These  low  values  of  R^  (compared 
to  the  critical  Rg  value  of  approximately  2300  for  pipe 
flow)  indicated  that  the  jet  flow  was  laminar.   The  nozzle 
was  then  temporarily  treated  as  a  circular  flow  channel. 
Kays  (36)  gives  the  length  to  diameter  ratio  necessary 
for  a  laminar  flow  velocity  profile  to  fully  develop  at  the 
entrance  of  a  circular  pipe  to  be 

length     Re 


diameter   20  (78) 

Thus,  for  the  case  under  consideration  these  ratios  had 
approximate  values  of  50,  which  was  many  times  larger 
than  the  nozzle  length  (measured  along  the  curvature  of 
the  nozzle)  to  diameter  ratio  of  2.4.   These  calculations 
were  interpreted  to  mean  that  the  flow  over  the  curved 
surface  of  the  nozzle  was  probably  not  fully  developed  by 
the  time  it  reached  the  exit.   Therefore,  the  boundary 
layer,  as  determined  experimentally,  should  exist. 

Schlichting  (_55)  gives  the  boundary  layer  thickness  6 
of  a  tube  with  infinite  radius  of  curvature  to  be 

^=-#1-,  (79) 

referred  to  the  tube  length  L^ .   The  length  of  the  nozzle 
was  2. 86  cm.   Based  on  this  model   ^  equals  2.95  mm  and 


2.S7  i;iin  for  Tests  I  and  II,  respectively.   Corrcspoiiding 
experimental  values  of  O     v/ere  2.5  mm  and  2.3  ;.i..., 
respectively.   Although  calculated  and  experimental  values 
differed  somewhat,  they  were  of  the  sane  order  oi      .     ...^ 
Thus,  the  experimental  boundary  ^„^,  „.  w.._^....^..  .^^^  .,„.-^ 
reasonable.   It  was  acknowledged  that  an  appropriate  value 
of  Re  was  dif:^'icul'c  to  obtain  since  the  velocity,  pressure, 
etc.,  were  changing  along  the  nozzle  length.   Hov/ever,  a 
more  sophisticated  theoretical  analysis  \vas  not  justified 
since  the  initial  condition  could  be  measured. 

Theoretical  models  for  velocity  profiles  near  nozzle 
exits,  such  as  those  of  Pai  (11)  and  Abramovicli  (12),  are 
usually  based  on  the  assumption  that  no  init        .dary 
layer  is  present  at  the  nozzle  exit.   This  p 
problem  in  trying  to  correlate  exijarimental  pro- 

files with  their  theoretical  -.odels .   To  correct  for  this 
situation  the  nondimensiw....^  ...._..:.  cocrdii..-„^  0^.0-^^.^.. ^ v. v- 
with  the  theoretical  model  was  modified  by  subtracting  a 
length  to  account  for  the  initial  bounaary  layer, 
magnitude  of  this  correction  was  Q.Q   mm  for  Test  I  and  5,8 
mm  for  Test  II.   Other  investigators  (11)-  (56)  have  made 

similar  corrections  for  initial  boundary  2:  ..■■^   ,. _,  _..., 

an  appropriate  length  to  the  experimental 

Comparison  of  theoretical  curves  wit. 
as  shown  in  Figures  11  and  16,  indicates  v 
velocity  decay  and  jet  spread  are  in  fairl 
with  that  proJ.icted  by  Pai's  (11)  theory. 

La.-.gmuii-  .  ..-.rements  .. 

the  nozzle  ex.  ^  .-..;.> ..^v..  ...o  elec 
essentially  uniform  across  the  jcv.  r.. 
not  only  for  the  pcontial  core  r^ 
mixing  a  ..  ,.   E.xperim> 

graphical^,        arcs  12  and  17. 


>  '..^i  fci  t«  ^4. -*  W    i^  <J 
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Upon  examination  of  the  experimental  data  two  conclu- 
sions become  immediately  clear.   First,  recombination  is 
essentially  negligible  within  the  jet  and  must  therefore 
occur  primarily  at  the  walls.   Secondly,  a  mechanism  must 
be  present  within  the  flow  which  accounts  for  the  fact 
that  the  electrons  along  a  given  radial  direction  possess 
a  constant  average  energy,  i.e.,  a  uniform  temperature. 
This  mechanism  can  best  be  understood  in  terms  of  the 
transport  processes  occurring  within  the  plasma  jet. 

As  has  been  pointed  out  previously,  the  electrons 
receive  most  of  the  field  energy  because  of  their 
relatively  small  mass  and  consequent  high  mobility.   This 
energy  is  transported  to  regions  remote  to  the  discharge 
zone  by  both  diffusion  and  conduction.   In  the  diffusion 
process  electrons  move  from  one  location  to  another  while 
transporting  finite  amounts  of  energy  to  the  new  location. 
During  the  diffusion  process  the  electrons  undergo  colli- 
sions with  other  particles  and  with  other  electrons.   These 
collisions  give  rise  to  energy  transport  by  conduction. 

The  electron-neutral  particle  collision  frequency  is 
very  high  compared  to  the  electron-electron  collision 
frequency.   However,  the  change  in  kinetic  energy  during  an 
electron-neutral  particle  collision  is  extremely  small  since 
a  more  effective  energy  transformation  occurs  between 
particles  of  equal  mass.   Thus,  the  energy  transferred  in 
electron-electron  collisions  becomes  important,  and  it  is 
this  mechanism  that  accounts  for  the  uniform  radial  electron 
temperatures.   The  importance  of  this  mechanism  was  also 
recognized  by  Gaither  (1_5)  in  his  investigation  of  the 
potential  core  region. 

In  Chapter  3  an  equation  was  derived  for  evaluating 
the  electron  number  density  as  a  function  of  r  and  z.   In 
order  to  make  use  of  this  equation  the  initial  electron 
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number  density  profile  f(r)  must  be  known,   ITor  pur-oscs 

of  this  investigation  it  was  observed  (see  Jigurc.^  _.  o.nd 
19)  that  f(r)  was  adequately  described  in  terms  of  a 
zeroth  order  Bessel  function  of  the  first  kind.   For  the 
best  curve  fit  of  the  initial  experimental  profiles  the 
factor   e~  1"^  ^■"'O  vvas  introduced  such  that 

^CTjs:^  *  ^o  ^"or   0:^ri=H«.       (80) 


■Rr;=0 


:or 


Note  that  Cj.  is  an  empirical  constant  determined  by  the 
initial  electron  number  density  profile.   C^  was  :     .ed  a 
value  of  141.  S  in  this  invest  ic-ation.   Th  ^  final  solution 
for  the  electron  nunibei- 


>-._.._>  —  ^        .,,->,       ,**V^i*    *-*»_- 


j2^ 

"'1-- 


where  n^^   is  the  electron  number  density  at  r  =  G      0. 
*-o 

Experimental  and  theoretical  profiles  for  axial  anv:: 
radial  electron  number  densities  ai-e  given  in  Figures  1-., 

15,  19,  and  20  for  tests  I  and  II.   Reason ., ,_^. . . 

is  observed.   Some  differences  between  theory  and  experi- 
ment were  expected  because  of  the  difficulxy  in  determinin.-; 
the  exact  breaking  point  on  the  Langmuir  probe  curves  and 
because  of  the  simplifications  introduced  in  deveL      the 
theory.   Certainly,  such  assumptions      .lowing  that  Do 
and  q^  were  constr^-.t  in  the  mixing      ..  wex-e  not 
desirable  but  wij.c  i.ucessary  in  or^..-    ■  ':  the  . 
not  becorae  extremely  complicated,   1\i^ 
such  assumptions  should  be  recog  .a-o 
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investigators  using  this  theory.   Further  discussion  of 
these  and  other  assumptions  are  given  in  Appendix  C. 

In  Test  III  a  potential  of  +71  volts  was  applied  to 
the  base  plate  in  the  base  I'egion.   This  particular  voltage 
level  was  chosen  because  field  distortion  effects  could  be 
introduced  at  this  level  while  holding  secondary  ioniza- 
tion effects  to  a  minimum.   The  experimental  results  are 
shown  graphically  in  Figures  21  through  24,   An  approximate 
increase  of  30%  over  Test  II  was  observed  in  the  electron 
temperature  at  the  nozzle  exit,  while  the  corresponding 
electron  number  density  was  decreased  by  appi^oximately  50%. 
The  superimposed  field  had  the  expected  effect  of  acceler- 
ating the  electrons,  resulting  in  a  higher  initial  electron 
temperature.   Also,  the  rate  of  electron  diffusion  was 
increased  such  that  a  lower  electron  number  density  was 
observed  along  the  jet  centerline.   Comparison  of  Figures 
19  and  23  shows  an  increase  in  the  electron  number  density 
in  the  base  region,  particularly  near  the  base  plate,  when 
a  potential  is  applied  to  thi§  region. 

It  was  concluded  from  this  test  (and  others  not  pre- 
sented herein)  that  the  electron  diffusion  process  in  a 
plasma  jet  flow  can  be  accelerated  or  retarded  by  applying 
a  potential  in  the  base  region.   The  extent  of  the  effect 
depends  on  the  potential  applied.   Large  gradients  are 
observed  in  the  plasma  potential  under  such  conditions, 
making  it  necessary  to  properly  bias  the  Langmuir  probes. 
Otherwise,  the  probes  will  be  damaged. 


6.   COXCLUSIONS  AND  RECOMMENDATIONS 

As  previously  stated, this  investigation  is  but  one 
step  toward  establishing'  a  cor.::lctc  solution  and  under- 
standing of  the  problem  ^^  ^ :.i...\_,v  v..  ^jo.rticio  bchavioy  ' 
low  pressure  plasma  jet  flows.   "    oxpeririiental  data 
presented  herein  have  been  shown  co  he   in  reasona',:      oc 
agreement  with  the  theory  employed  and  with        ^rimental 
results  of  other  investigators.   Tr.ase  data  support  the 
following  conclusions  for  the  type  of  flov/  inv.     -tad: 

1.  Radial  electron  temperature  distribution^  >--vs 
constant  Throughout  the  potential  core,  jet 
mixing,  and  base  regions. 

2.  Very  little  recombination  occurs  within  the  initial 
region  of  the  jet.   Tr.iLj  ii  verified  by  the  fact 
that  the  electron  temp-^^ ..  ^..res  v^...^^..   -..._.  ._/^.__ 
high,  particularly  in  'che  radial  directiov. . 

3.  Radial  electron  number  density  profii.      rhe 
initial  region  of  the  jet  are  described  quite 

by  the  theory  presented  herein,  provided  that  thv^ 
initial  profile  is  properly  chosen. 

4.  The  need  for  taking  into  account      „.  w..  .,^^...., 
i^nc    thermal  conduction  is  clearly  _ 

5.  The  electron  number  density  and  cl. 

profiles  are  continuous  in  the  i      _  boundary 
layer  at  the  nozzle  exit.   The  electron  tompcraturo 
was  obiorved  to  remain  constant  throu^      this 
lay.  .• . 

6.  An  ...>.. x-ud  potential  in  t-ne  uus^. 

.;lectrical  characteristics  : 

in  the-  jet, 
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This  writer  has  found  the  plasma  jet  to  be  an  extremely 
valuable  research  tool  and  recommends  the  following  areas 
of  research  for  future  investigations: 

1.  The  extension  of  this  investigation  to  the 
supersonic  regime. 

2.  The  extension  of  this  study  to  flows  having  an 
initial  turbulent  mixing  region. 

3.  The  devising  of  a  series  of  experiments,  similar 
in  nature  to  the  one  considered  herein,  with  the 
objective  of  determining  the  effects  of  pressure 
over  a  wide  range  (low  micron  range  to  atmospheric) 
on  electron  gas  behavior, 

4.  The  development  of  techniques,  particularly 
electronic,  which  would  permit  rapid  analysis  of 
Langmuir  probe  characteristics. 

5.  Further  examination  of  the  effects  of  applied 
fields  on  transport  properties. 
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FIGURK  1 .   PLASMA  GENERATOR 
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FIGURE    5.      PLASMA   GEIIERATOR   AIID   TEST    SECTION 
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FIGURE  6.   PLASMA  JET 
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FIGURE  12.   RADIAL  ELECTRON  TEMPERATURE  DISTRIBUTION 
IN  ARGON  PLASMA  JET  FOR  TEST  I 
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FIGURE  14.   RADIAL  ELECTRON  DENSITY  DISTRIBUTION 
IN  ARGON  PLASMA  JET  FOR  TEST  I 
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FIGURE  17.   R.^DIAL  ELECTRON  TEMPER-ATURE  DISTRIBUTION 
IN  ARGON  PLASMA  JET  FOR  TEST  II 
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RADIAL  ELECTRON  DENSITY  DISTRIBUTION 
IN  ARGON  PLASMA  JET  FOR  TEST  II 
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FIGURE  21.   RADIAL  ELECTRON  TEMPERATURE  DISTRIBUTION 
IN  ARGON  PLASMA  JET  FOR  TEST  III 
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APPENDICES 


A.       SUMMARY   OF   PliYSi^ii.   D-.Tii    i:Oxi  ^niuo.^    ^^i) 


Symbol  itoni 


M     Molecular  weight 


nij^  Aiass   of   atom 


Value 


3S.944 


Kff 


mole 
6.63    X    10-26   Kg 


Ref  eroncc- 


I'j^  Radius   of   atom 


1.32  X  10-10  meters  (53) 


R 


Gas  constant 


20S.1 


Joules 


»v-i— A 


(5S) 


^p    Specific  heat  at 
constant  pressure 


523 


(5S) 


'y    Specific  hoai;  a' 
constant  volume 


;15  - 


(ni:) 


k     Ratio  of  s,.;... 

heats 


Sutherland  constant 


K     Coefficient  of  tnermal   0.0162  — "  ~"'       "" < 

'        J   ^  •  .  me'ter-scc-A 

conductivity 


^-s)- 


Vi    First  ionization 
potential 


15.3  ev 


(27)  (2^ 


Vj^^    Second  ionization 
potential 


27.  _ 


^Civcn   for  T   -   -73   K  and  ? 


^   atiaospnero 
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B.   THEORETICAL  CONSIDERATIONS  FOR  TKE 
DETERiAIINATION  OF  MEAN  FREE  PATHS 


In  the  plasma  flow  system  uxider  consideration  each 
xype  of  particle  collides  wit--      ether  types  of 
particles  as  well  as  with  its  ov/n  type.   Thus,  since  those 
interactions  have  a  bearing  on  probe  characteristics  and 
continuum  considerations,  it  is  of  interest  to  know  the 
mean  free  paths  involved.   A  method   for  determining  these 
mean  free  paths  has  been  outlined  by  Graf  (13)  and  is 
given  in  the  following  paragraphs. 

For  neutral-neutral  particle  interactions  the  /.;oan 
free  path  is  given  by  Sears  (5  7)  as 


where  Vj^  is  the  viscosity,  T  is  the  temperature,  and  P  is 
the  pressure.   French  (51)  gives  a  graphical         .\:c.- 
tion  of  this  equation  based  on  the  Ler.nard-Jones  6-12 

potential  model  for  the  calcula". „...  „.... .. ..y. 

Assuming  the  momentum  cross  sect-on  for  neutral-ion 
encounters  to  oe  the  same  as  for  neutral-neutral  encounter 
the  cross  sections  Q  are  then  given  by 


The  collisions  considered 


V  J.  ■«—     --  W  -. 


transfer   with   a  mean   doflcction   of   90   dercreos . 

o 

"^It   must   bo    ..    .hat   this    _..   ....„   ..-  ^., 

since  charge  oxchauj-'O  i.iay  occur.   However,        oxchn 

cross  section  data  are  rather  scarce  at  1 

elation  of  e.xisting  data  suggest 

of  neutral-^on  encounters  are  of 

as  tho.-,r:  for  momentum  transfer  (.  .  .  . 
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ss 


^-v  =  ^^-^  ^ '^-a;  •  CB3) 

Substituting  for  Qi_j  in  equation  (16)  gives 


(85) 


Thus, 

The   neutral-ion  r;ican   free  path   is   given  by- 


TV..  .=   -^  ._    -  _.  (--) 


Hence, 


The  electron-ion  mean  free  path  is  given  by  ti:e 
equation 

A  _•  =  — '- — ^  \c-^j 

where  v^  is  the  mean  electron  velocity,  and      __•  is 
the  collision  frequency  for  momentum  transfer, 
expression  for   lK_--    ^*^s  been  derived  by  Dreice:.'  (5S) 
and  is  given  by 


^rrc^^-/^/- 


^  -  ■  ■ ' ,' 


where  n^   is  the  number  of  particles  in  a  sphere  of  radius 

equal  to  one  Dcbye  length,*^  and   £^   is  the  perraittivity 

of  free  space. 

From  expressions  for  the  various  relaxation  times 

between  the  plasma  components,  Graf  (IG)  obtains  the 

following  expressions  for  the  mean  free  paths  of  the  charged 

comoonents  in  terms  of  A   . : 

c-i 


^c-e=  ^.-:.  '--^la  (90) 


.,.e  =  A,_,  (rcjr^i^   •  ^ 


(SI) 


A,_;  =  A  _■  .  (92) 


For  electron-neutral  particle  encounters  ...o  ...^an  fre; 
path  is  determiiied  from  -co  equation 


?\   ,  ~  r .  (93  ) 


'.^v-v^^^v_/A.    \.-l.^'<b  -..v 


A  graphical  represo...c^on  of  Q^.^-^  fo:.  ^^ .. 
(less  tnan  1.8  ev)  is  given  in  Figuiv  ' 

No  expression  was  given  by  Graf  ,  ^_^ . 

Ay^_Q   should  be  larger  than   /^e-n  ^*''      °^   ^^^^   electron 
density  being  lower  than  the  neutral  particle  density. 


,  Debye  length  is  that  c 
plasma  system  beyond  which  th. 
test  particle  exert  negligible  ^... 


C.   ESTIMATE  OF  UNCERTAINTIES  INVOLVED  IN  THE 
MEASUREMENT  OF  ELECTRON  NUMBER  DENSITIES 


It  is  necessary  in  an  investigation  of  this  type  to 
establish  whenever  possible  the  extent  of  the  uncertainties 
introduced  due  to  assumptions  and  methods  of  measurement. 
Such  an  analysis  then  forms  a  basis  for  a  comparison 
between  experimental  and  theoretical  results. 

Wehner  and  Medicus  (59)  have  written  an  excellent 
paper  on  possible  sources  of  error  involved  in  making 
probe  measurements.   The  most  important  sources  of  error 
mentioned  in  their  paper  which  have  a  direct  bearing  on 
the  work  done  in  this  investigation  may  be  enumerated  as 
follows : 

1.  Disturbance  of  plasma  by  probe. 

2.  Probe  contamination. 

3.  Electron  emission  by  probe  due  to  photon,  ion, 
and  metastable  atom  impact. 

4.  Non-Maxwellian  distribution  of  electrons. 

5.  Inability  to  measure  plasma  potential  accurately 
(and  therefore  electron  density) . 

The  first  source  of  error  was  minimized  by  making  the 
Langmuir  probe  as  small  as  possible.   Thus,  the  electron 
diffusion  losses  to  the  probe  were  minimized.   The  degree 
of  disturbance  of  the  plasma  by  the  probe  was  checked  by 
observing  the  effect  of  the  measuring  probe  on  a  neigh- 
boring probe.   At  a  distance  of  0.5  mm  no  appreciable 
effect  was  observed.   Thus,  the  disturbance  due  to  the 
probe  was  considered  to  be  negligible. 
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9  J. 

Probe  cor.'..^~.:.^-.c^^:^v, ..  ^...^...r^ ^^.\. 

are  "two  of  tlio  mosv.    '  ,o  ovcrccnio  ' .. 

making  probe  measure:n>-  a 

given  probe  chan        :.c  shov/c  stic 

changed  as  a  function  of  time.   Cleaning  t^   .  .  .U 

the  original  characteristic.   Thu~ .  it  %vaF;  nc 

clean  the  probo  \^^,.^.>j^    .tv-^-w^'^^.-..  ^. .i  ^^rw.-'v:;  \,.\:^-^^ .      v,v^..- 

taniination  eifeccs  were  then  con£;idered  tc  .aaii. 

Hov/ever,  little  could  be  done  about  the      aary  emission 
effects.   For  a  low  density  argon  plas'^:.  .16) 

estimated  the  experirr.ent 

mately  4%  too  high  due  to  j...^^^    ■ ^ ..o..-.   o 

amount  of  error  soc/.c-v,  ..o  ije   a  reasonable  estimate  for 
the  case  under  consideration. 

Graphical  plots  ox  ln(JQ    )    versus  Vv-,  for  the  retai^ding 
field  region  of  the  probe  characteristics  gave  iresulv 
straight  lines.   As  was  previously  explained  in  section  2,6, 

the  resulting  straigh'^  '^ .  a.ndic:  .^v..  „....„  .      ..^^^-  ... 

distribuxicn  was  inaeed  Maxwellian.   Thus,  no  ^.r^^^le:.:  was 
encountered  in  this  investigation  due  to  non-Ma;vWeliian 
distributions , 

The  problem  in       ..vostiga'.     .hich  ^.     ■   ont 

the  greatest  difficulty  w  :S  in  detc _,  exactly  ■-- 

the  "break"  occurred  in  v.-^j  iOiiii-log  plot  of  'che  pi-obe 
characteristic.   The  conventional  "two-tangent  ii^itorccpt" 
method,  as  shown  in  Figure  10,  was  used  herein,  and  the 
intercept  of  the  tv/o  tangents  was  used  to  locate  the 
"break."   The  diff icv.ljy  v/as  that  the  curves  did  not  bruak 
sharply,  and  above  ■^' '-•  so-. --""or,  ^.\.. 

non-linear.   Graf  ('.  .,      v.,^  ^^„u^  „^.. .;  >,a..  .v-^cs 

may  give  results  that 
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As  was  previously  mentioned,  the  drift  in  the  electron 
current  density  was  sometimes  found  to  be  as  large  as  10% 
due  to  the  long  recording  times.   This  variation  should 
have  occurred  mainly  between  the  measurements  of  the 
electron  density  profiles  rather  than  within  a  given 
profile,  i.e.,  the  time  between  the  measurement  of  profiles 
was  usually  longer  than  that  required  for  measuring  a  given 
profile . 

Thus,  the  combined  effects  of  the  above  sources  of 
error  allow  for  errors  as  great  as  30%  in  the  measurement 
of  the  electron  numbei-  densities.   Generally  speaking,  the 
measured  value  will  be  too  high  because  of  the  method  used 
in  interpreting  the  probe  characteristics. 

The  reproducibility  of  the  electron  temperature  was 
determined  to  be  within  about  15  to  20%  of  a  given  level  for 
a  given  set  of  conditions.   Variations  in  electron  tempera- 
ture due  to  making  semi-log  plots  from  photographs  of 
oscilloscope  traces  was  estimated  to  be  approximately  300 
to  500  K,  regardless  of  the  temperature  level.   Such 
differences  show  up  in  the  determination  of  the  electron 
number  density  as  the  inverse  of  the  square  root  of  the 
electron  temperature. 

In  making  comparisons  between  the  measured  values  and 
the  predicted  electron  density  values,  the  reader  should 
also  be  aware  that  some  very  basic  assumptions  have  been 
made  in  arriving  at  a  theoretical  model.   For  example,  D 
and  q   were  assumed  to  be  constant  in  the  region  of  interest, 
i.e.,  near  the  nozzle  exit.   The  effects  of  these  assump- 
tions will  be  magnified  near  the  edge  of  the  plasma  flow 
and  at  locations  far  removed  from  the  nozzle.   Perhaps  a 

more  realistic  approach  in  future  investigations  would  be 

N 
go  try  to  use  mean  values  for  D^^  and  q^  or  to  take  their 

variation  into  account. 
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